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1. Introduction 



The large-A limit of Gauge Theory was proposed in 1974 by 't Hooft [1] 
for quantum chromo dynamics (QCD). The dimensionality of the gauge 
I group SU(Ac) was used as a parameter, considering Ac as a large num- 

ber and performing an expansion in 1/Nc- The motivation was an expan- 
sion in the inverse number of field-components A in statistical mechanics 
where it is known as the 1/A-expansion, and is a standard method for 
non-perturbative investigations. 

The expansion of QCD in 1/Ac rearranges diagrams of perturbation 
theory in a way which is consistent with a string picture of strong inter- 
action, whose phenomenological consequences agree with experiment. The 
accuracy of the leading-order term, which is often called large-A QCD or 
multicolor QCD, is expected to be of the order of the ratios of meson widths 
to their masses, i.e. about 10-15%. 
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While QCD is simplified in the large- A^c limit, it is still not yet solved. 
Generically, it is a problem of infinite matrices, rather than of infinite vec- 
tors as in the theory of second-order phase transitions in statistical me- 
chanics. 

We start these Lectures by showing how the 1/A'^-expansion works for 
the 0(Ar)-vector models, and describing some applications to the four-Fermi 
interaction and the nonlinear sigma model. Then we concentrate on the 
Yang-Mills theory at large Nc- 

The methods described in these Lectures are developed mostly in the 
seventies and the beginning of the eighties. They are used over and over 
again when discussing a relation between Gauge Theory and String Theory. 

The conformal invariance, described in Section 2 for the four-Fermi 
interaction, is used in the modern approaches based on the AdS/CFT cor- 
respondence [2] (discussed in the lectures by L. Thorlacius at this School). 
The notion of the large-iV limit, described in Section 3, is widely used, in 
particular, in Matrix Theory [3] (discussed in the lectures by W. Taylor 
at this School). The loop equation, described in Section 4, has been re- 
cently applied [4] for studying the AdS/CFT correspondence. The large- 
reduced models, described in Section 5, are applied to a nonperturbative 
formulation of superstring [5] and very recently to noncommutative gauge 
theory [6]. 

The content of my fifth lecture at this School, which was devoted to an 
application of the large- iV methods to Matrix Theory at finite temperature, 
is not included in the text. It is mostly contained in Ref. [7]. 

2. 0{N) Vector Models 

The simplest models, which become solvable in the limit of a large number 

of field components, deal with a field which has N components forming 
an 0{N) vector in an internal symmetry space. A model of this kind was 
first considered by Stanley [8] in statistical mechanics and is known as 
the spherical model. The extension to quantum field theory was done by 
Wilson [9] both for the four-Fermi and 99^ theories. 

In the framework of perturbation theory, the four-Fermi interaction 
is renormalizable only in d = 2 dimensions and is non-renormalizable 
for d > 2. The 1/iV-expansion resums perturbation-theory diagrams after 
which the four-Fermi interaction becomes renormalizable to each order in 
for 2 < d < 4. An analogous expansion exists for the nonlinear 0(N) 
sigma model. 

The 1/iV expansion of the vector models is associated with a resumma- 
tion of Feynman diagrams. A very simple class of diagrams — the bubble 
graphs — survives to the leading order in 1/N. This is why the large- A?" 
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limit of the vector models is solvable. Alternatively, the large- A*" solution 
is nothing but a saddle-point solution in the path-integral approach. The 
existence of the saddle point is due to the fact that A'' is large. This is to be 
distinguished from a perturbation-theory saddle point which is due to the 
fact that the coupling constant is small. Taking into account fluctuations 
around the saddle-point results in the 1/A^-cxpansion of the vector models. 

We begin this Section with a description of the 1/A^-expansion of the 
A^-component four-Fermi theory analyzing the bubble graphs. Then we in- 
troduce functional methods and construct the l/A/"-expansion of the 0{N)- 
symmetric nonlinear sigma model. At the end we discuss the factorization 
in the 0{N) vector models at large A^. 

2.1. FOUR-FERMI INTERACTION 

The action of the 0(A^)-symmetric four-Fermi interaction in a d- 
dimensional Euclidean space is defined by 

S[iP,iP]= J d'^x(^^dip + mipip-^{iijij;fy (1) 

Here d = '^^d^ and ^ = (fpi, ■ ■ ■ ^^Pn) is a spinor field which forms an 
A'-component vector in an internal symmetry space so that 

N 

i^i; = Y,^i^l;i. (2) 
1=1 

In d = 2 this model was studied in the large-AT limit in Ref. [10] and is 
often called the Gross-Neveu model. 

The dimension of the four-Fermi coupling constant G is dim [G] = 
j^2-d Pqj. ^]^^g j.gg^gon, the perturbation theory for the four-Fcrmi inter- 
action is renormalizable in d = 2 but is non-renormalizable for d > 2 (and, 
in particular, in ci = 4). This is why the old Fermi theory of weak interac- 
tions was replaced by the modern electroweak theory, where the interaction 
is mediated by the and Z bosons. 

The action (1) can be equivalently rewritten as 

S[^,^,X] = J d''x(^i;d^ + m^^-x^^ + ^y (3) 

where x is an auxiliary field. The two forms of the action, (1) and (3), are 
equivalent due to the equation of motion 



X = Gil^iJ: 



(4) 
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which can be derived by varying the action (3) with respect to x- 

In the path-integral quantization, where the partition function is defined 

by 

Z = j DxDi)Di)e-^\^^^^^\ (5) 

with S ^, x] given by Eq. (3), the action (1) appears after performing 
the Gaussian integral over x- Therefore, one alternatively gets 

Z = j Di^Dil^e-^i^^-^] (6) 

with S [V', V'] given by Eq. (1). 

The perturbative expansion of the 0(iV)-symmetric four- Fermi theory 
can be conveniently represented using the formulation (3) via the auxiliary 
field X- Then the diagrams arc of the type of those in Yukawa theory, and 
resemble the ones for QED with and being an analog of the electron- 
positron field and x being an analog of the photon field. However, the 
auxiliary field x(a^) does not propagate, since it follows from the action (3) 
that 

i^o {x-y)^ {X{x)x{y) )Gauss = G5^<'\x - y) (7) 

or 

Ab) = (x(-p)x(p))Gauss = G (8) 
in momentum space. 

It is convenient to represent the four-Fermi vertex as the sum of two 
terms 




(9) 



where the empty space inside the vertex is associated with the propaga- 
tor (7) (or (8) in momentum space). The relative minus sign makes the 
vertex antisymmetric in both incoming and outgoing fermions as is pre- 
scribed by the Fermi statistics. 

The diagrams that contribute to second order in G for the four-Fermi 
vertex are depicted, in these notations, in Figure 1. The O(A^) indices 
propagate through the solid lines so that the closed line in the diagram in 
Figure lb corresponds to the sum over the 0(iV) indices which results in a 
factor of N . Analogous one- loop diagrams for the propagator of the '0-field 
are depicted in Figure 2. 

Remark on the one-loop Gell-Mann-Low function of four- Fermi theory 
Evaluating the diagrams in Figure 1 which are logarithmically divergent 
in d = 2, and noting that the diagrams in Figure 2 do not contribute to 
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d) e) 

Figure 1. Diagrams of the second order of perturbation theory for the four- Fermi vertex. 
The diagram b) involves the sum over the 0(N) indices. 




a) b) 

Figure 2. One-loop diagrams for the propagator of the ■i/'-field. The diagram b) involves 
the sum over the 0{N) indices. 

the wave- function renormalization of the ■0- field, which emerges to the next 
order in G, one gets for the one-loop Gell-Mann-Low function 

si,a)^-^Ji^. (10) 

The four-Fcrmi theory in 2 dimensions is asymptotically free as was first 
noted by Anselm [11] and rediscovered in Ref. [10]. 

The vanishing of the one-loop Gell-Mann-Low function in the Gross- 
Neveu model for N = lis related to the same phenomenon in the Thirring 
model. The latter model is associated with the vector-like interaction 
(V'7/iV') of species of fermions with 7^ being the 7-matrices in 2 di- 
mensions. Since a bispinor has in d = 2 only two components V'l and ^^2, 
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Figure 3. Bubble diagram which survives the large- A/" limit of the 0{N) vector models. 

both the vector-hke and the scalar-hke interaction (1) for A?^ = 1 reduce 
to ?/'iV'i'02V'2 since the square of a Grassmann variable vanishes. Therefore, 
these two models coincide. For the Thirring model, the vanishing of the 
G ell- Mann-Low function for any G was shown by Johnson [12] to all loops. 



2.2. BUBBLE GRAPHS AS ZEROTH ORDER IN 1/N 

The perturbation-theory expansion of the 0(A/^)-symmetric four-Fermi the- 
ory contains, in particular, the diagrams of the type depicted in Figure 3 
which are called bubble graphs. Since each bubble has a factor of N, the 
contribution of the n-bubble graph is oc G^+^AT" which is of the order of 

G^+^AT" ~ G (11) 

as iV — DO since 

G-^. (12) 

Therefore, all the bubble graphs are essential to the leading order in 1/N. 
Let us denote 

n loops 

= G+...+ G^(^+ G^+^(^---(^ +... 

(13) 

In fact the wavy line is nothing but the propagator D of the x fi^ld with 
the bubble corrections included. The first term G on the RHS of Eq. (13) 
is nothing but the free propagator (8). 

Summing the geometric series of the fermion-loop chains on the RHS of 
Eq. (13), one gets analytically^ 

1 r rl'^h sp [k + im] [k + f) + im\ 

^Recall that the free Euclidean fermionic propagator is given by So(p) = (ip + m)~^ 
due to Eqs. (3), (5) and the additional minus sign is associated with the fermion loop. 
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Figure 4- Some diagrams of the l/iV-expansion for the 0{N) four- Fermi theory. The 
wavy line represents the (infinite) sum of the bubble graphs (13). 



This determines the exact propagator of the x field at large N. It is O {N~^) 
since the couphng G is included in the definition of the propagator. 

The idea is now to change the order of summation of diagrams of per- 
turbation theory using 1/N rather than G as the expansion parameter. 
Therefore, the zeroth-order propagator of the expansion in 1 /N is defined 
as the sum over the bubble graphs (13) which is given by Eq. (14). Some of 
the diagrams of the new expansion for the four-Fermi vertex axe depicted in 
Figure 4. The first diagram is proportional to G while the second and third 
ones are proportional to G^ or G^, respectively, and therefore are of order 
O(iV-i) or 0(iV-2) with respect to the first dia gram. The perturbation 
theory is thus rearranged as the l/A?^-expansion. 

The general structure of the 1/A^-expansion is the same for all vector 
models, say, for the A'^-component nonlinear sigma model which is consid- 
ered in Subsection 2.4. 

The main advantage of the expansion in 1/A?^ for the four- Fermi inter- 
action, over the perturbation theory, is that it is renormalizable in d < 4 
while the perturbation-theory expansion in G is renormalizable only in 
d = 2. Moreover, the 1/A^-expansion of the four- Fermi theory in 2 < d < 4 
demonstrates [9] an existence of an ultraviolet-stable fixed point, i.e. a 
nontrivial zero of the Gell-Mann-Low function. 

In order to show that the 1/A^-expansion of the four- Fermi theory is 
renormalizable in 2 < c? < 4, let us analyze indices of the diagrams of the 
l/AT-expansion. First of all, we shall get rid of an ultraviolet divergence of 
the integral over the d-momentum k in Eq. (14). The divergent part of the 
integral is proportional to A'^^'^ (logarithmically divergent in d = 2) with 
A being an ultraviolet cutoff. It can be canceled by choosing 



N 



(15) 
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where is a proper dimensionless constant which is not necessarily positive 
since the four- Fermi theory is stable with either sign of G. The power of A in 
Eq. (15) is consistent with the dimension of G. This prescription works for 
2 < d < 4 where there is only one divergent term while another divergency 
oc p^lnA emerges additionally in d = 4. This is why the consideration is 
not applicable in c? = 4. 

The propagator D{p) is therefore finite, and behaves at large momenta 
\p\ 2> m as 

D{P) ^ ^ • (16) 

The standard power-counting arguments then show that the only divergent 
diagrams appear in the propagators of the ■)/' and x fields, and in the ■i/'-x-V' 
three-vertex. These divergencies can be removed by a renormalization of 
the coupling g, mass, and wave functions of and x- 

This completes a demonstration of renormalizability of the 1/A^-expan- 
sion for the four-Fermi interaction in 2 < d < 4. For more detail, see 
Ref. [13]. 

2.3. SCALE AND CONFORMAL INVARIANCE OF FOUR-FERMI THEORY 

The coefficient in Eq. (15) can easily be calculated in d = 3. To evaluate 
the divergent part of the integral in Eq. (14), we put p = and m = 0. 
Remembering that the 7-matrices are 2x2 matrices in d = 3, we get 

d^k spkk _^ d?k 1 _ 1 _ A 

J (2^ w"" y j2^¥~^J "'i'"!-^- 

Note that the integral is linearly divergent in d = 3 and A is the cutoff 
for the integration over \k\. This divergence can be canceled by choosing G 
according to Eq. (15) with g equal to 

g*=7r. (18) 

To calculate in d = 3 the coefficient of proportionality in Eq. (16), let us 
choose G = tt'^/NA as is prescribed by Eqs. (15), (18) and put in Eq. (14) 
m = since we are interested in the asymptotics \p\ S> m. Then the RHS 
of Eq. (14) can be rearranged as 

The integral is obviously convergent. 

Equation (19) (or (16) in d dimensions) is remarkable since it shows 
that the scale dimension of the field % changes its value from = d/2 in 



/c^ -I- kp 
k'^{k+p)'^ 



1 



N\p\ 



(19) 
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a) b) 

Figure 5. Diagrams for the 1/iV-correction to the ip-S.e\d. propagator (a) and the 
three- vertex (b). 



perturbation theory to Z,^ = 1 in the zeroth order of the 1/N expansion 
(remember that the momentum-space propagator of a field with the scale 
dimension I is proportional to This appearance of scale invariance 

in the l/A?^-expansion of the four-Fermi theory at 2 < d < 4 was first 
pointed out by Wilson [9] and implies that the Gell-Mann-Low function 
B{g) has a zero at g = g^.. which is given in d = 3 by Eq. (18). 

The (logarithmic) anomalous dimensions of the fields ip, Xj and of the 
V'-X-f/' three- vertex in d = 3 to order 1/N can be found as follows. The 
1/iV-correction to the propagator of the ^j-ficld is given by the diagram 
depicted in Figure 5a). Since we are interested in an ultraviolet behavior, 
we can put again m = 0. Analytically, we have 



^ d?k k+p 
(27r)3 \k\{k+pf ■ 



(20) 



The (logarithmically) divergent contribution emerges from the domain of 
integration \k\ 3> \p\ so we can expand the integrand in p. The p-indepen- 
dent term vanishes after integration over the directions of k so that we 
get 



S ^(j)) = ip 



1 + 



8 d-^k 1 
N J (27r)3 



ip 



1 + 



37r2Ar 



In + finite 



(21) 

The diagram, which gives a non-vanishing contribution to the three- 
vertex in order 1/A^, is depicted in Figure 5b. It reads analytically 



8 d^k {k + Pi){k+P2) 



i2nf\k\ik+pi)\k + p2)^ 



(22) 



where pi and p2 the incoming and outgoing fermion momenta, respectively. 
The logarithmic domain is ^ |p|max with \p\niax being the largest of \pi\ 
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and \p2\- This gives 

2 A2 

r(pi , P2) = 1 - In ^— + finite . (23) 

An analogous calculation of the correction for the field X is a bit 
more complicated since involves three two-loop diagrams (see Ref. [14]). 
The resulting expression for D~^{p) reads 



{nd{p)) ' = 



A 
9^ 



1 

+ 



/2 A2 ^ 
2A - IpI ( - In + finite 



' (24) 

The linear divergence is canceled to order \jN providing g is equal to 

5* = TT f 1 + 1 V (25) 



iV. 

which determines to order 1/iV. After this 0"^^^) takes the form 

To make all three expressions (21), (23), and (26) finite, we need loga- 
rithmic renormalizations of the wave functions of ■)/'- and x-fields and of the 
vertex F. This can be achieved by multiplying them by the renormalization 
constants ^ 

Z,(A) = l-7iln^ (27) 

where /i stands for a reference mass scale and 7^ arc anomalous dimensions. 
The index i stands for ^, or v for the ■i/'- and x-propagators or the three- 
vertex r, respectively. We have, therefore, calculated 



7V 



37r2Ar' 
2 



to order l/AT. Due to Eq. (4) 7-^ coincides with the anomalous dimension 
of the composite fields •^i^) 

li,^ = Ix ■ (29) 

Note, that 



ZlZ-^Z^ = l. (30) 
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This implies that the effective charge is not renormahzed and is given by 
Eq. (25). Thus, the nontrivial zero of the Gell-Mann-Low function persists 
to order \/N (and, in fact, to all orders of the l/A'^-expansion). 

If g is chosen exactly at the critical point 5* , then the renormalization- 
group equations 

d In r j 



7i 



{9') , (31) 



where Tj stands generically either for vertices or for inverse propagators, 
possess the scale invariant solutions 

Vi oc . (32) 

For the four-Fermi theory in cZ = 3, Eq. (32) yields 

Sip) = hi^'" ^ (33) 



' 2 \ 
p ^ 



rto.P.) = KTfKni^^]. (35) 
\PiJ \Pi Pi J 

where / is an arbitrary function of the dimensionless ratios which is not 
determined by scale invariance. The indices here obey the relation 

7v = 7^- + ^7x (36) 

which guarantees that Eq. (30), implied by scale invariance, is satisfied. 

The indices 7, are given to order 1/N by Eqs. (28). When expanded in 
1/A^, Eqs. (33) and (34) obviously reproduces Eqs. (21) and (26). Therefore, 
one gets the exponentiation of the logarithms which emerge in the 1/N- 
expansion. The calculation of the next terms of the 1/iV-expansion for the 
indices 7.^ is contained in Rcf. [15]. 

Scale invariance implies, in a renormalizable quantum field theory, more 
general conformal invariance as is first pointed out in Refs. [16, 17]. The 
conformal group in a d-dimensional space-time has {d+l){d + 2)/2 param- 
eters as is illustrated by Table 1. More about the conformal group can be 
found in the lecture by Jackiw [18]. 

A heuristic proof [16] of the fact that scale invariance implies conformal 
invariance is based on the explicit form of the conformal current K^, which 
is associated with the special conformal transformation, via the energy- 
momentum tensor: 

= (2x,x" - x^6'^) . (37) 
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Group 



Transformations 



# Parameters 



Lorcntz 


rotations 




d(d-l) 
2 


Poincare 


+ d translations 




d{d+l) 
2 


Wcyl 


+ 1 dilatation 


= p x^ 


d^+d+2 
2 


Conformal 


+ d special conformal 




(d+l)(d+2) 
2 



TABLE 1. Contents and the number of parameters of groups of space-time 
symmetry. 



Differentiating, we get 



(38) 



which is proportional to divergence of the dilatation current. Therefore, 
both the dilatation and conformal currents vanish simultaneously when 
9fj_i, is traceless which is provided, in turn, by the vanishing of the Gell- 
Mann-Low function. 

Conformal invariance completely fixes three-vertices as was first shown 
by Polyakov [19] for scalar theories. The proper formula for the four-Fermi 
theory reads [20] 



r( ^ 2>£(l)£(|_7v) 
d'^k k + pi 



k +P2 



(39) 



where the coefficient in the form of the ratio of the F-functions is prescribed 
by the normalization (33) and (34) and the indices are related by Eq. (36) 
but can be arbitrary otherwise^. 



^The only restriction 7y, > is imposed by the Kallen-Lehmann representation of the 
propagator while there is no such restriction on 7^ since it is a composite field. 
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Equation (39), which results from conformal invariance, unambiguously 
fixes the function / in Eq. (35). In contrast to infinite-dimensional confor- 
mal symmetry in d = 2, the conformal group in d > 2 is less restrictive. It 
fixes only the tree-point vertex while, say, the four-point vertex remains an 
unknown function of two variables. 

The integral on the RHS of Eq. (39) looks in d = 3 very much like 
that in Eq. (22) and can easily be calculated to the leading order in 1/N 
when only the region of integration over large momenta with > |p|max = 
max{|pi|, \p2\} is essential to this accuracy. 

Let us first note that the coefficient in front of the integral is oc 7v ~ 
1/N, so that one has to peak up the term ~ l/j^ in the integral for the 
vertex to be of order 1. This term comes from the region of integration with 
1^1 k, bimax- Recalling that \pi - p2\ < \p\max in Euclidean space, one gets 

f (fk k+pi k+p2 1 1 



27r [{k + pi)2]l+7x/2 [(fc + p2)2]l+7x/2 \J,\l+2^^-^j2 (^2^^) 



While the integral in Eq. (40) is divergent in the ultraviolet for 7v < 0, this 
divergence disappears after the renormalization. 

Equation (23) is reproduced by Eq. (40) when expanding in 1/N. This 
dependence of the three-vertex solely on the largest momentum is typical 
for logarithmic theories in the ultraviolet region where one can put, say, 
Pi = without changing the integral with logarithmic accuracy. This is 
valid if the integral is fast convergent in infrared regions which is our case. 

Remark on broken scale invariance 

Scale (and conformal) invariance at a fixed point g = g* holds only for 
large momenta \p\ 3> m. For smaller values of momenta, scale invariance 
is broken by masses. In fact, any dimensional parameter /x breaks scale 
invariance. If the bare coupling g is chosen in the vicinity of 5*, then scale 
invariance holds even in the massless case only for |p| 2> /x where g{p) 
approaches g^ while it is broken if \p\ < fi. 




(40) 



where Eq. (36) has been used and 




(41) 
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2.4. NONLINEAR SIGMA MODEL 

The nonlinear 0{N) sigma model^ in 2 Euclidean dimensions is defined by 
the partition function 

where n = {rii, . . . , n^v) is an O(A^) vector. While the action in Eq. (42) is 
pure Gaussian, the model is not free due to the constraint 

n\x) = ^, (43) 

which is imposed on the n field via the (functional) delta-function. 

The sigma model in d = 2 is sometimes considered as a toy model for 
QCD since it possesses: 

1) asymptotic freedom [21]; 

2) instantons for AT = 3 [22]. 

The action in Eq. (42) is ~ A'' as A'^ ^ oo but the entropy, i.e. a contribu- 
tion from the measure of integration, is also ^ N so that a straightforward 
saddle point is not applicable. 

To overcome this difficulty, we introduce an auxiliary field u{x), which 
is ~ 1 as AT — DO, and rewrite the partition function (42) as 



Z cx J Du{x) J Dn{x) e \- V 



(44) 



where the contour of integration over u(x) is parallel to imaginary axis. 
Doing the Gaussian integration over n, we get 

Z (X Du{x) e 2 ^ ^'J2^J ^ ' , (45) 

The first term in the exponent is nothing but the sum of one-loop diagrams 
in 2 dimensions 



^Svln(-dl + u{x)) = . M . (46) 

^The name comes from elementary particle physics where a nonlinear sigma model in 
4 dimensions is used as an effective Lagrangian for describing low-energy scattering of 
the Goldstone 7r-mesons. 
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where the auxihary field u is denoted again by the wavy fine. 

Now the path integral over u{x) in Eq. (45) is a typical saddle-point 
one: the action ~ N while the entropy ~ 1 since only one integration over 
u is left. The saddle-point equation for the nonlinear sigma model 

\-NG{x,x-u^^) = {) (47) 



9 

while G is defined by 



G{x,y-u) = ( y 



-dl + u 



(48) 



The coupling in Eq. (47) is ~ as is prescribed by the con- 

straint (43) which involves a sum over N terms on the LHS. This guaran- 
tees that a solution to Eq. (47) exists. Next orders of the 1/iV-expansion 
for the 2-dimensional sigma model can be constructed analogously to the 
previous Subsections. 

The l/A'^-expansion of the 2-dimensional nonlinear sigma model has 
many advantages over perturbation theory, which is usually constructed 
solving explicitly the constraint (43), say, choosing 



1 

riN = -A 



N-l 

l-g^Y^nl (49) 

a=l 



and expanding the square root in . Only iV — 1 dynamical degrees of free- 
dom are left so that the 0(A'^)-symmetry is broken in perturbation theory 
down to 0{N — 1). The particles in perturbation theory are massless (like 
Goldstone bosons) and it suffers from infrared divergencies. 
On the contrary, the solution to Eq. (47) has the form 

nsp = m| = A2e"^, (50) 

where A is an ultraviolet cutoff. Therefore, all A'^ particles acquires the same 
mass mR in the l/A^-expansion so that the 0{N) symmetry is restored. This 
appearance of mass is due to dimensional transmutation which says in this 
case that the parameter mR, rather than the renormalized coupling constant 
g'^ is observable. The emergence of the mass cures the infrared problem. 

To show that (50) is a solution to Eq. (47), let us look for a transla- 
tionally invariant solution Usp(x) = Tn\. Then Eq. (47) in the momentum 
space reads 

1 (Pk 1 N f^' dk"^ AT A2 



g"^ J (27r)2 k"^ + m^^ 47r Jo k"^ + m\ Att m^^ 
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The exponentiation results in Eq. (50). 

Equation (51) relates the bare coupling and the cutoff A and allows 
us to calculate the Gell-Mann-Low function yielding 

The analogous one- loop perturbation-theory formula for any N reads [21] 

(53) 

Thus, the sigma-model is asymptotically free in 2-dimensions for A?^ > 2 
which is the origin of the dimensional transmutation. There is no asymp- 
totic freedom ioi N = 2 since 0(2) is Abelian. 



2.5. LARGE-A/' FACTORIZATION IN VECTOR MODELS 

The fact that a path integral has a saddle point at large N implies a very 
important feature of large-A'^ theories — the factorization. It is a general 
property of the large- A^ limit and holds not only for the O(A^) vector mod- 
els. However, it is useful to illustrate it by these solvable examples. 

The factorization at large N holds for averages of singlet operators, for 
example 

( u (xi) ...u (xfe) ) = Z-^ j Due -f Spi-[-a2H+3^ /<^^-" ^ (^^) . . . ^ (^^) 

^ (54) 
in the 2-dimensional sigma model. 

Since the path integral has a saddle point at some u [x) = Usp (x) which 
is, in fact, a;-independent due to translational invariance, we get to the 
leading order in 1/N: 

{ u (xi) ...u {xk) ) = Usp (xi) ...Usp (xjfc) -I- O (^-^^ , (55) 

which can be written in the f actor ized form 

{uixi)...u{xk)) = {u{xi))...{u{xk))+0 1^^^ . (56) 

Therefore, u becomes "classical" as A^ — > oo in the sense of the 1/N- 
expansion. This is an analog of the WKB-expansion in h = 1/N. "Quan- 
tum" corrections arc suppressed as 1/N. 

We shall return to discussing the large- A^ factorization in the next Sec- 
tion when considering the large- A^ limit of QCD. 
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3. Large-N QCD 

The method of the 1/iV-expansion can be apphed to QCD. This was done by 
't Hooft [1] using the inverse number of colors for the gauge group SU(A'^c) 
as an expansion parameter. 

For a SU(A''c) gauge theory without virtual quark loops, the expansion 
goes m l/Nj and rearranges diagrams of perturbation theory according 
to their topology. The leading order in l/AT^ is given by planar diagrams, 
which have a topology of a sphere, while the expansion in 1/iV^ plays the 
role of a topological expansion. This reminds an expansion in the string 
coupling constant in string models of the strong interaction, which also has 
a topological character. 

Virtual quark loops can be easily incorporated in the l/A/'j^-expansion. 
One distinguishes between the 't Hooft limit when the number of quark 
flavors Nf is fixed as iVc oo and the Veneziano limit [23] when the 
ratio Nf/Nc is fixed as iVc ^ oo. Virtual quark loops arc suppressed in 
the 't Hooft limit as 1/Nc and lead in the Veneziano limit to the same 
topological expansion as dual-resonance models of strong interaction. 

The simplification of QCD in the large- AT^ limit is due to the fact that 
the number of planar graphs grows with the number of vertices only expo- 
nentially rather than factorially as do the total number of graphs. Correla- 
tors of gauge invariant operators factorize in the large- A^c limit which looks 
like the leading-order term of a "semiclassical" WKB-expansion in l/Nc- 

We begin this Section with a description of the double-line represen- 
tation of diagrams of QCD perturbation theory and rearrange it as the 
topological expansion in l/Nc- Then we discuss some properties of the 
l/ATc-expansion for a generic matrix- valued field. 

3.1. INDEX OR RIBBON GRAPHS 

In order to describe the l/A'c-expansion of the Yang-Mills theory, it is 
convenient to represent the gauge field by a Hermitean N x N matrix 



Here [t"]*-^ are the generators of the gauge group (a = 1, . . . — 1 for 
SU (iVc)) with the normalization 




(57) 



a 



2 



(58) 



The (Euclidean) action reads 




(59) 
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where 

F^^u = d^A^, - duAf, - i[Af„ A^] (60) 

is the non-Abehan field strength and g is the couphng constant. 
The propagator of the matrix field A^^{x) has the form 

4 (X) At' (,) = i (SH^' - If'S"') (.-,), (61) 

where we have assumed, as usual, a gauge fixing to define the propagator 
in perturbation theory. For instance, one has 

in the Feynman gauge. 

Equation (61) can be immediately derived from the standard formula 

A1{x)Al{y)) =5'^'D,,{x-y) (63) 

/ Gauss 

multiplying by the generators of the SU(A''c) gauge group according to the 
definition (57) and using the completeness condition 

1 1 



^ {ey^ {f'f = -t: ^^^^^ - —S'^S''^ for SU(iVe) , (64) 



a=l ^ ^ 

where the factor of 1/2 is due to the normalization (58) of the generators. 

Wc concentrate in this Section only on the structure of diagrams in the 
index space, i.e. the space of the indices associated with the SU(iVc) group. 
We shall not consider, in most cases, space-time structures of diagrams 
which are prescribed by Feynman's rules. 

Omitting at large Nc the second term in parentheses on the RHS of 
Eq. (61), we depict the propagator by the double line 

A^^ (x) Al^ (y) ) oc g'^S'^S''^ = ^. — T" f ■ (65) 

Each line represents the Kronecker delta-symbol and has orientation which 
is indicated by arrows. This notation is obviously consistent with the space- 
time structure of the propagator which describes a propagation from x to 

y- 

The arrows are due to the fact that the matrix is Hermitean and 
its off-diagonal components are complex conjugate. The independent fields 
are, say, the complex fields for i > j and the diagonal real fields AJ^. 
The arrow represents the direction of the propagation of the indices of the 
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complex field vl^ for i > j while the complex-conjugate one, A^^ = (vl^)*, 
propagates in the opposite direction. For the real fields A^^, the arrows are 
not essential. 

The double-line notation appears generically in all models describing 
matrix fields in contrast to vector (in internal symmetry space) fields whose 
propagators are depicted by single lines as in the previous Section. 

The three- gluon vertex, which is generated by the action (59), is de- 
picted in the double-line notations as 




(66) 

where the subscripts 1, 2 or 3 refer to each of the three gluons. The rel- 
ative minus sign is due to the commutator in the cubic in A term in the 
action (59). The color part of the three- vertex is antisymmetric under in- 
terchanging the gluons. The space-time structure, which reads in the mo- 
mentum space as 

= <^M1M2 iPi - P2)f,s + <^/X2M3 iP2 - P3)^i + (^Mi/^3 (P3 - Pl)^2 , (67) 

is antisymmetric as well. We consider all three gluons as incoming so that 
their momenta obey ^1+^2+^3 = 0. The full vertex is symmetric as is 
prescribed by Bose statistics. 

The color structure in Eq. (66) can alternatively be obtained by multi- 
plying the standard vertex 

bl>P2,P3) = /"^"^"^7^iM2^3 (Pl,P2,P3) (68) 

by (t"i)^'^' (t"2)^2J2 (t«3)»3i3^ ^j^j^ fabc i^gj^g ^Yi_e structure constants of the 
SU(iVc) group, and using the formula 

jaxa2a3 ^^aiyiji fj.a2y2j2 ^^asysjs _ 1 ^(Jiii3(ji2ii^i3i2 _ (Jiii2(ji3ii(ji2i3 j ^ 

^ (69) 
which is a consequence of the completeness condition (64). 

The four-gluon vertex involves six terms — each of them is depicted 
by a cross — which differ by interchanging of the color indices. We depict 
the color structure of the four-gluon vertex for simplicity in the case when 
h = J2 = i, 12 = jz = i, h = ji = k, = ji = / but i,j,k,l take on 
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Figure 6. Double-line representation of a one-loop diagram for the gluon propagator. The 
sum over the A^c indices is associated with the closed index line. The relative contribution 
of this diagram is ~ g^N^ ~ 1. 



different values. Then only the following term is left 



TP 



(70) 



and there are no deltas on the RHS since the color structure is fixed. In other 
words, we pick up only one color structure by equaling indices pairwise. 

The diagrams of perturbation theory can now be completely rewritten in 
the double-line notation [1]. The simplest one which describes the one-loop 
correction to the gluon propagator is depicted in Figure 6.^ This diagram 
involves two three-gluon vertices and a sum over the Nc indices which is 
associated with the closed index line. Therefore, the relative contribution 
of this diagram is ~ g^Nc- 

In order for the large- A'^c limit to be nontrivial, the bare coupling con- 
stant should satisfy 

5^ ~ — • (71) 

This dependence on Nc is similar to Eqs. (12) and (47) for the vector models 
and is prescribed by the asymptotic-freedom formula 

2 _ 247r^ 
^ " 117V,ln(A/AQCD) ^ ^ 

of the pure SU(A^c) gauge theory. 

Thus, the relative contribution of the diagram of Figure 6 is of order 

Figure 6 ~ g^Nc ~ 1 (73) 

in the large- A^c limit. 



*Here and in the most figures below the arrows of the index lines are omitted for 
simplicity. 
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The double lines of the diagram in Figure 6 can be viewed as bounding 
a piece of a plane. Therefore, these lines represent a two-dimensional object 
rather than a one-dimensional one as the single lines do in vector models. 
These double-line graphs are often called in mathematics the ribbon graphs 
or fatgraphs. We shall see below their connection with Riemann surfaces. 

Remark on the U{Nc) gauge group 

As is said above, the second term in the parentheses on the RHS of Eq. (64) 
can be omitted at large Nc- Such a completeness condition emerges for the 
V{Nc) group whose generators {A = 1, . . . , N^) are 

= 73^) ' t.T^T^ = \5^^. (74) 
They obey the completeness condition 

A=l ^ 

The point is that elements of both the SU(A'^2) group and the \]{Nc) group 
can be represented in the form U = expiB, where i? is a general Hermitean 
matrix for \J{Nc) and a traceless Hermitean matrix for SU(A''c). 

Therefore, the double-line representation of perturbation-theory dia- 
grams which is described in this Section holds, strictly speaking, only for 
the U(A^c) gauge group. However, the large-A^c limit of both the U(A^c) 
group and the S'U{Nc) group is the same. 

3.2. PLANAR AND NON-PLANAR GRAPHS 

The double-line representation of perturbation theory diagrams in the in- 
dex space is very convenient to estimate their orders in l/Nc- Each gluon 
propagator contributes a factor of g'^ Each three- or four-gluon vertex con- 
tributes a factor of g^'^. Each closed index line contributes a factor of Nc- 
The order of g in 1/Nc is given by Eq. (71). 

Let us consider a typical diagram for the gluon propagator depicted 
in Figure 7. It has eight three-gluon vertices and four closed index lines 
which coincides with the number of loops. Therefore, the relative order of 
this diagram in 1/Nc is 

Figure 7- (fif^jYc)^ ~ 1. (76) 

The diagrams of the type in Figure 7, which can be drawn on a sheet 
of a paper without crossing any lines, are called the planar diagrams. For 



for U(iVe) 



(75) 



22 



YURI MAKEENKO 




Figure 7. Double-line representation of a four-loop diagram for the gluon propagator. 
The sum over the A^c indices is associated with each of the four closed index lines whose 
number is equal to the number of loops. The contribution of this diagram is ~ g^N* ~ 1. 




Figure 8. Doublo-liiic representation of a three-loop non-planar diagram for the gluon 
propagator. The diagram has six three-gluon vertices but only one closed index line (while 
three loops!). The order of this diagram is ~ g^Nc ~ 

such diagrams, an adding of a loop inevitably results in adding of two three- 
gluon (or one four-gluon) vertex. A planar diagram with n2 loops has n2 
closed index lines. It is of order 



n2-loop planar diagram 




so that all planar diagrams survive in the large- A'^c limit. 

Let us now consider a non-planar diagram of the type depicted in Fig- 
ure 8. This diagram is a three-loop one and has six three-gluon vertices. 
The crossing of the two lines in the middle does not correspond to a four- 
gluon vertex and is merely due to the fact that the diagram cannot be 
drawn on a sheet of a paper without crossing the lines. The diagram has 
only one closed index line. The relative order of this diagram in l/Nc is 

Figure 8 ~ g^N^ ~ . (78) 

It is therefore suppressed at large Nc by 

The non-planar diagram in Figure 8 can be drawn without line-crossing 
on a surface with one handle which is usually called in mathematics a torus 
or the surface of genus one. A plane is then equivalent to a sphere and 
has genus zero. Adding a handle to a surface produces a hole according 
to mathematical terminology. A general Riemann surface with h holes has 
genus h. 
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The above evaluations of the order of the diagrams in Figures 6-8 can 
now be described by the unique formula 



Thus, the expansion in 1/Nc rearranges perturbation-theory diagrams ac- 
cording to their topology [1]. For this reason, it is referred to as the topo- 
logical expansion or the genus expansion. The general proof of Eq. (79) for 
an arbitrary diagram is given in Subsection 3.3. 

Only planar diagrams, which are associated with genus zero, survive 
in the large- AT^ limit. This class of diagrams is an analog of the bubble 
graphs in the vector models. However, the problem of summing the planar 
graphs is much more complicated than that of summing the bubble graphs. 
Nevertheless, it is simpler than the problem of summing all the graphs, 
since the number of the planar graphs with no vertices grows at large no 
exponentially [24, 25] 



while the number of all the graphs grows with no factorially. There is no 
dependence in Eq. (80) on the number of external lines of a planar graph 
which is assumed to be much less than uq. 

It is instructive to see the difference between the planar diagrams and, 
for instance, the ladder diagrams which describe e^e~ elastic scattering in 
QED. Let the ladder has n rungs. Then there are n! ladder diagrams, but 
only one of them is planar. This simple example shows why the number 
of planar graphs is much smaller than the number of all graphs, most of 
which are non-planar. 

We shall discuss in these Lectures what is known about solving the 
problem of summing the planar graphs. 

Equation (79) holds, strictly speaking, only for the relative order while 
the contribution of tree diagrams to a connected n-point Green's function is 
~ {g'^)"'~^ which is its natural order in l/Nc- In order to make contributions 
of all planar diagrams to be of the same order ~ 1 in the large- A^c limit, 
independently of the number of external lines, it is convenient to contract 
the Kronecker deltas associated with external lines. 

Let us do this in a cyclic order as is depicted in Figure 9 for a generic 
connected diagram with three external gluon lines. The extra deltas which 
are added to contract the color indices are depicted by the single lines. They 
can be viewed as a boundary of the given diagram. The actual size of the 
boundary is not essential — it can be shrunk to a point. Then a bounded 
piece of a plane will be topologically equivalent to a sphere with a puncture. 




genus 



(79) 



#p (no) = # of planar graphs ~ const' 



no 



(80) 
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Figure 9. Generic index diagram with no = 10 vertices, n\ — 10 gluon propagators, 
712 = 4 closed index lines, and B = 1 boundary. The color indices of the external lines 
are contracted by the Kronecker deltas (represented by the single lines) in a cyclic order. 
The extra factor of 1/Nc is due to the normalization (81). Its order in 1/Nc is ~ l/iV^ 
in accord with Eq. (79). 



I shall prefer to draw planar diagrams in a plane with an extended boundary 
(boundaries) rather than in a sphere with a puncture (punctures). 

It is clear from the graphic representation that the diagram in Figure 9 
is associated with the trace over the color indices of the three-point Green's 
function 

G^ms ^2, X3) = ]^ ( tr (^l) {X2) A^, (X3)]) . (81) 

We have introduced here the factor 1/Nc to make of 0{1) in the large- 
Nc limit. Therefore, the contribution of the diagram in Figure 9 having one 
boundary should be divided by Nc- 

The extension of Eq. (81) to multi-point Green's functions is obvious: 

ixi,...,Xn)^^JtT [A^, (xi) ■■■A,^ (x„)]) . (82) 

The factor 1/Nc, which normalizes the trace, provides the natural normal- 
ization 

= 1 (83) 

of the averages. 

Though the two terms in the index-space representation (66) of the 
three-gluon vertex look very similar, their fate in the topological expansion 
is quite different. When the color indices are contracted anti-clockwise, the 
first term leads to the planar contributions to G^^^ , the simplest of which is 
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a) b) 

Figure 10. Planar a) and non-planar b) contributions of the two color structures in 
Eq. (66) for three-gluon vertex to G^^^ in the lowest order of perturbation theory. 



depicted in Figure 10a. The anti-clockwise contraction of the color indices 
in the second term leads to a non-planar graph in Figure 10b which can be 
drawn without crossing of lines only on a torus. Therefore, the two color 
structures of the three-gluon vertex contribute to different orders of the 
topological expansion. The same is true for the four-gluon vertex. 

Remark on oriented Riemann surfaces 

Each line of an index graph of the type depicted in Figure 9 is oriented. 
This orientation continues along a closed index line while the pairs of index 
lines of each double line have opposite orientations. The overall orientation 
of the lines is prescribed by the orientation of the external boundary which 
we choose to be, say, anti-clockwise like in Figure 10a. Since the lines are 
oriented, the faces of the Riemann surface associated with a given graph 
are oriented too — all in the same way — anti-clockwise. Vice versa, such 
an orientation of the Riemann surfaces unambiguously fixes the orientation 
of all the index lines. This is the reason why we omit the arrows associated 
with the orientation of the index lines: their directions are obvious. 

Remark on cyclic- ordered Green's functions 

The cyclic-ordered Green's functions (82) naturally arise in the expansion of 
the trace of the path-ordered non-Abelian phase factor for a closed contour. 
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V 





a) 



b) 



Figure 11. Example of a) connected and b) disconnected planar graphs. 



One gets 



( 



1 



) 



n=0 p j.^ 



/ 



Hi— tin 



'in) 



(xi, . . . , Xfi) . 



(84) 



The reason is because the ordering along a closed path implies the cyclic- 
ordering in the index space. 

Remark on generating Junctionals for planar graphs 

By connected or disconnected planar graphs we mean, respectively, the 
graphs which were connected or disconnected before the contraction of the 
color indices as is illustrated by Figure 11. The graph in Figure 11a is 
connected planar while the graph in Figure lib is disconnected planar. 

The usual exponential relation between the generating functionals W [J\ 
and Z [J] for connected graphs and all graphs, does not hold for the planar 
graphs. The reason is that an exponentiation of such a connected planar 
diagram for the cyclic-ordered Green's functions (82) can give disconnected 
non-planar ones. 

The generating functionals for all and connected planar graphs can 
be constructed [26] by means of introducing non-commutative sources 
j^{x). "Non-commutative" means that there is no way to transform 
j'/ii (^i)-^'m2 (^2) into 3 ^2^x2)3 ^^{xi). This non-commutativity of the sources 
reflects the cyclic-ordered structure of the Green's functions (82) which 
possess only cyclic symmetry. 

Using the short-hand notations 




(85) 
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Figure 12. Graphic derivation of Eq. (89): Z [j] is denoted by an empty box, W [j] is 
denoted by a shadow box, j is denoted by a filled circle. By picking a leftmost external line 
of a planar graph, we end up with a connected planar graph, whose remaining external 
lines are somewhere to the right interspersed by disconnected planar graphs. It is evident 
that jZ [j\ plays the role of a new source for the connected planar graph. If we instead 
pick up the rightmost external line, we get the inverse order Z[j]j, which results in 
Eq. (90). 



where the symbol o includes the sum over the d- vector (or whatever avail- 
able) indices except for the color ones, we write down the definitions of the 
generating functionals for all planar and connected planar graphs, respec- 
tively, as 

Z[j]^f:(^tH3oAr) (86) 

n=0 ^^^'^ I 

and 

^[j] = E(i^tr(j°^r) ■ (87) 

n=Q \ c / conn 

The planar contribution to the Green's functions (82) and their con- 
nected counterparts can be obtained, respectively, from the generating func- 
tionals Z [j] and W [j] applying the non-commutative derivative which is 
defined by 

TT^j.(y) / (j) = V^^'^ - y) / (i) > (88) 

where / is an arbitrary function of j's. In other words the derivative picks 
up only the leftmost variable. 

The relation which replaces the usual one for planar graphs is 

Z\j\ = W\^Z\ji\\, (89) 

while the cyclic symmetry says 

W\jZ\j\\ = W{Z\ji\j\. (90) 

A graphic derivation of Eqs. (89) and (90) is given in Figure 12. In other 
words, given W [j], one should construct an inverse function as the solution 
to the equation 

j^{x) = J^{x)W\j\, (91) 
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after which Eq. (89) says 

Z[j] = W[J]. (92) 

More about this approach to the generating functionals for planar graphs 
can be found in Ref. [27]. 

An iterative sohition to Eq. (91) for the Gaussian case can be easily 
found. In the Gaussian case, only G^^^ is nonvanishing which yields 

W[j] = l-g''joDoj, (93) 

where the propagator D is given by Eq. (62). Using Eq. (89), we get ex- 
plicitly 

Z [j] = 1-9^ j d^xd^yD^^^x - y) j ^{x) Z [j] 3,{y) Z [j] . (94) 

While this equation for Z [j] is quadratic, its solution can be written only 
as a continued fraction due to the non-commutative nature of the variables. 
In order to find it, we rewrite Eq. (94) as 

^^^"^^ l + g^!ddx ddy D,4x - y) j^{x) Z [j] j,{y) ' ^^^^ 
whose iterative solution reads [26] 

^ = ^ ^^7^^ • (96) 

1 + ^^' 1 ^ 

^ + 9 3 -^D^3 

^ + 9 3 3 



3.3. TOPOLOGICAL EXPANSION AND QUARK LOOPS 

It is easy to incorporate quarks in the topological expansion. A quark field 
belongs to the fundamental representation of the gauge group SU(A^c) and 
its propagator is represented by a single line 

(X 5ij = i — ' — j . (97) 

The arrow indicates, as usual, the direction of propagation of a (complex) 
field ip. We shall omit these arrows for simplicity. 

The diagram for the gluon propagator which involves one quark loop is 
depicted in Figure 13a. It has two three gluon vertices and no closed index 
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a) b) 

Figure 13. Diagrams for the gluon propagator with a quark loop which is represented 
by the single lines. The diagram a) involves one quark loop and has no closed index lines 
so that its order is ~ ~ '^/N^- The diagram b) involves three loops one of which is a 
quark loop. Its relative order is ~ g^N"^ ~ 1/iVc. 



lines so that its order in 1/Nc is 

Figure 13a ~ 5^ ~ . (98) 



N. 



c 



Analogously, the relative order of a more complicated tree-loop diagram in 
Figure 13b, which involves one quark loop and two closed index lines, is 



Figure 13b ~ g^N^ ~ . (99) 



c 



It is evident from this consideration that quark loops are not accom- 
panied by closed index lines. One should add a closed index line for each 
quark loop in order for a given diagram with L quark loops to have the 
same double-line representation as for pure gluon diagrams. Therefore, 
given Eq. (79), diagrams with L quark loops are suppressed at large Nc 

by 

(j^ \ L+2-genus 

The single-line representation of the quark loops is similar to the one of 

the external boundary in Figure 9. Moreover, such a diagram emerges when 
one calculates perturbative gluon corrections to the vacuum expectation 
value of the quark operator 

Or = ^^rV;, (101) 
where T stands for one of the combinations of the gamma-matrices: 

r = I, 75, 7m, ^^75, S^i/ = ^.[lt^,lu] ■ (102) 
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The factor of 1/N(. is introduced in (101) to make it 0{1) in the large- 
Nc hmit. Therefore, the external boundary can be viewed as a single line 
associated with valence quarks. The difference between virtual quark loops 
and external boundaries is that each of the latter comes along with the 
factor of l/Nc due to the definitions (82) and (101). 

In order to prove Eqs. (79) and its quark counterpart (100), let us 

(3) 

consider a generic diagram in the index space which has Hq three-point 
vertices (either three-gluon or quark-gluon ones), four-gluon vertices, 
ni propagators (either gluon or quark ones), n2 closed index lines, L virtual 
quark loops and B external boundaries. Its order in 1/Nc is 

^_(^2-jm-4^'-n^'*'jY"2 ^ ]^n2-ni+no-B {'^^'^) 

where the total number of vertices no = nQ^-* + is introduced. The extra 
factor of is due to the extra normalization factor of l/N^ in operators 

associated with external boundaries. 

The exponent on the RHS of Eq. (103) can be expressed via the Euler 
characteristics x of a given graph of genus h. Let us first mention that a 
proper Riemann surface, which is associated with a given graph, is open 
and has B + L boundaries (represented by single lines). This surface can be 
closed by attaching a cap to each boundary. The single lines then become 
double lines together with the lines of the boundary of each cap. We have 
already considered this procedure when deducing Eq. (100) from Eq. (79). 

The number of faces for a closed Riemann surface constructed in such 
a manner is n2 + L + B, while the number of edges and vertices are ni and 
no, respectively. Euler's theorem says that 

X = 2-2h = n2 + L + B-ni+no. (104) 

Therefore the RHS of Eq. (103) can be rewritten as 

jY"2-ni+no-S _ j^2-2h-L-2B (105) 

We have thus proven that the order in 1 /N^ of a generic graph does not 
depend on its order in the coupling constant and is completely expressed 
via the genus h and the number of virtual quark loops L and external 
boundaries B by 

/ -j^ X 2/i+L+2(B-l) 

generic graph ~ yiV~J ' (-'-^^) 



For S = 1, we recover Eqs. (79) and (100). 
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Remark on the order of gauge action 

We see from Eq. (82) that the natural variables for the large- A^c limit are the 
matrices which inchidc the factor of g (see Eq. (57)). In these variables, 
the action (59) is 0{N^) at large N^, since g'^ is ~ 1/Nf. and the trace is 

This result can be anticipated from the free theory because the kinetic 
part of the action involves the sum over — 1 free gluons. Therefore, the 
non-Abelian field strength (60) is ~ 1 for g^ ~ l/Nc- 

The fact that the action is O(N^) in the large- A^c limit is a generic prop- 
erty of the models describing matrix fields. It will be crucial for developing 
saddle-point approaches at large Nc which are considered below. 

Remark on phenomenology of multicolor QCD 

While Nc = Sin the real world, there are phenomenological indications that 
1/Nc may be considered as a small parameter. We have already mentioned 
some of them in the text — the simplest one is that the ratio of the p- 
meson width to its mass, which is ~ l/Nc, is small. Considering l/Nc 
as a small parameter immediately leads to qualitative phenomenological 
consequences which are preserved by the planar diagrams associated with 
multicolor QCD, but are violated by the non-planar diagrams. 

The most important consequence is the relation of the l/A?^c-expansion 
to the topological expansion in the dual-resonance model of hadrons. Vast 
properties of hadrons are explained by the dual-resonance model. A very 
clear physical picture behind this model is that hadrons are excitations of 
a string with quarks at the ends. 

I shall briefly list some consequences of multicolor QCD: 

1) The "naive" quark model of hadrons emerges at A'^c = oo. Hadrons 
are built out of (valence or constituent) quark and antiquark qq, while 
exotic states like qqqq do not appear. 

2) The partial width of decay of the ^-meson, which is built out of ss 
(the strange quark and antiquark), into K~^K~ is ~ l/A^o while that 
into 7r+7r~7r'' is ~ l/A'c- This explains Zweig's rule. The masses of the 
p- and w-mesons are degenerate at A^c = oo. 

3) The coupling constant of meson-meson interaction is small at large Nc- 

4) The widths of glucballs are ~ l/A^^ , i-^- they should be even narrower 
than mesons built out of quarks. The glueballs do not interact or mix 
with mesons i^i Nc = oo. 

All these hadron properties (except the last one) approximately agree 
with experiment, and were well-known even before 1974 when multicolor 
QCD was introduced. Glueballs are not yet detected experimentally (maybe 
because of their property listed in the item 4). 
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3.4. LARGE-A^c FACTORIZATION 

The vacuum expectation values of several colorless or white operators, 
which are singlets with respect to the gauge group, factorize in the large- 
Nc limit of QCD (or other matrix models). This property is similar to that 
already discussed in Subsection 2.5 for the vector models. 

The simplest gauge-invariant operator in a pure SU(A^c) gauge theory 
is the square of the non-Abelian field strength: 

0(x) = ^trF2,(x). (107) 

The normalizing factor is the same as in Eqs. (81), (82), which provides the 
natural normalization 




(108) 



The contribution of all planar graphs to the average on the LHS of Eq. (108) 
is of order 1 in accord with the general formula (106) for B = 1. 

In order to verify the factorization in the large- A^c limit, let us consider 
the index space diagrams for the average of the product of two colorless op- 
erators O (xi) and O (^2) given by (107). It involves a factorized part when 
gluons are emitted and absorbed by the same operators. The contribution 
of the factorized part is or order 1 as above. 

Alternatively, the connected correlator of the two operators is associated 
with the general formula (106) for two boundaries B = 2. Its contribution 
is suppressed by l/N^ in the laige-Nc limit. For this correlator, at least 
one gluon line is emitted and absorbed by different operators O (xi) and 
O (X2). Notice, that these graphs themselves are planar, while the suppres- 
sion comes from the number of boundaries. 

This example illustrates the general property that only (planar) dia- 
grams with gluon lines emitted and absorbed by the same operators survive 
as Nc ^ 00. Since correlations between the colorless operators O (,ti) and 
O (X2) are of order the factorization property holds as ^ 00: 

For a general set of gauge-invariant operators Oi, . . . , On, the factor- 
ization property can be represented by 

(Oi---On) = (Oi)---(0„)+o(^) . (110) 
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This is analogous to Eq. (56) for the vector models. 

The factorization in large-iVc QCD was first discovered by A. A. Migdal 
in the late seventies. An important observation that the factorization im- 
plies a semiclassical nature of the large- A^c limit of QCD was done by Wit- 
ten [28]. We shall discuss this in the next two Subsections. 

The factorization property also holds for gauge-invariant operators con- 
structed from quarks like in Eq. (101) as a consequence of Eq. (106). 

Remark on factorization beyond perturbation theory 

The large- A^c factorization (110) has been shown above to all orders of per- 
turbation theory. It can be also verified at all orders of the strong coupling 
expansion in the SU(A^c) lattice gauge theory. A non-perturbative proof 
of the factorization will be given in the next Section by using quantum 
equations of motion (the loop equations). 

3.5. THE MASTER FIELD 

The large-A^c factorization in QCD assumes that gauge-invariant objects 
behave as c-numbers, rather than as operators. Likewise the vector models, 
this suggests that the path integral is dominated by a saddle point. 

We already saw in Subsection 2.5 that the factorization in the vector 
models does not mean that the fundamental field itself, for instance ft in 
the sigma-model, becomes "classical". It is the case, instead, for a singlet 
composite field. 

We are now going to apply a similar idea to the Yang-Mills theory 
whose partition function reads 



The action, ~ N^, is large as Nc ^ oo, but the "entropy" is also ~ due 
to the — 1 integrations over A": 



Consequently, the saddle-point equation of the large- A^c Yang-Mills theory 
is not the classical one. 

The idea is to rewrite the path integral over for the Yang-Mills 
theory as that over a colorless composite field <I> [A], likewise it was done 
in Subsection 2.4 for the sigma-model. The expected new path-integral 
representation of the partition function (111) would be something like 




(111) 




(112) 




(113) 
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The Jacobian 



in Eq. (113) is related to the old entropy factor, so that J[$] ~ 1 in the 
large- A^c limit. 

The original partition function (111) can be then rewritten as 

Z(xyi^$e^cM*]-iv|5[*], (115) 

where S [$] represents the Yang-Mills action in the new variables. The new 
"entropy" factor is 0{1) because the variable $ [A\ is a color singlet. 
The large parameter enters Eq. (115) only in the exponent. Therefore, 
the saddle-point equation can be immediately written: 

Remembering that <^ is a functional of A^: $ = $ [^] , we rewrite the 
saddle-point equation (116) as 

'& = (V„F,.r= (117) 



It differs from the classical Yang-Mills equation by the term on the RHS 
coming from the Jacobian (114). 

Given J [$] which depends on the precise from of the variable $ [A] , 
Eq. (117) has a solution 

A,{x) = A';lix). (118) 

Let us first assume that there exists only one solution to Eq. (117). Then 
the path integral is saturated by a single configuration (118), so that the 
vacuum expectation values of gauge-invariant operators are given by their 
values at this configuration: 

{0) = 0{a^^{x)). (119) 

The factorization property (110) will obviously be satisfied. 

An existence of such a classical field configTiration in multicolor QCD 
was conjectured by Witten [28]. It was discussed in the lectures by Cole- 
man [29] who called it the master field. Equation (117) which determines 
the master field is often referred to as the master-field equation. 

A subtle point with the master field is that a solution to Eq. (117) is de- 
termined only up to a gauge transformation. To preserve gauge invariance, 
it is more reasonable to speak about the whole gauge orbit as a solution of 
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Eq. (117). However, this will not change Eq. (119) since the operator O is 
gauge invariant. 

The conjecture about an existence of the master field has surprisingly 
rich consequences. Since vacuum expectation values are Poincare invari- 
ant, the RHS of Eq. (119) does. This implies that A'^^{x) must itself be 
Poincare invariant up to a gauge transformation: a change of A^^(a^) under 
translations or rotations can be compensated by a gauge transformation. 
Moreover, there must exist a gauge in which ^^'(x) is space-time indepen- 
dent: ^^'(x) = ^^^(0) . In this gauge, rotations must be equivalent to a 
global gauge transformation, so that ^^^(0) transforms as a Lorentz vector. 

In fact, the idea about such a master field in multicolor QCD may 
be incorrect as was pointed out by Haan [30]. The conjecture about an 
existence of only one solution to the master-field equation (117) seems to 
too strong. If several solutions exists, one needs an additional averaging 
over these solutions. This is a very delicate matter, since this additional 
averaging must still preserve the factorization property. One might better 
think about this situation as if ^^'(0) would be an operator in some Hilbert 
space rather than a c-valued function. Such an operator-valued master field 
is sometimes called the master field in the weak sense, while the above 
conjecture about a single classical configuration of the gauge field, which 
saturates the path integral, is called the master field in the strong sense. 

The concept of the master field is rather vague until a precise form of 
the composite field <I? [A], and consequently the Jacobian $ [A\ that enters 
Eq. (117), is not defined. However, what is important is that the master 
field (in the weak sense) is space-time independent. This looks like a simpli- 
fication of the problem of solving large- A^c QCD. A Hilbert space, in which 
the operator A^^{Q) acts, should be specified by <5 [A\. We shall consider in 
the next Subsection a realization of these ideas for the case of $ [A\ given 
by the trace of the non-Abelian phase factor for closed contours. 

Remark on non-commutative probability theory 

An adequate mathematical language for describing the master field in mul- 
ticolor QCD (and, gencrically, in matrix models at large N^) was found by 
I. Singer in 1994. It is based on the concept of free random variables of 
non-commutative probability theory, introduced by Voiculescu [31]. How 
to describe the master field in this language and some other applications of 
non-commutative free random variables to the problems of planar quantum 
field theory are discussed in Refs. [32, 33]. 
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3.6. l/Nc AS SEMICLASSICAL EXPANSION 

A natural candidate for the composite operator $ [A\ from the previous 
Subsection is given by the trace of the non-Abehan phase factor for closed 
contovirs — the Wilson loop. It is labeled by the loop C in the same sense 
as the field {x) is labeled by the point x, so we shall use the notation 

$(C) = $ [A] = i-trPe'i'c'^^''"^''^^). (120) 

Nobody up to now managed to reformulate QCD at finite in terms 
of <5 (C) in the language of path integral. This is due to the fact that self- 
intersecting loops are not independent (they are related by the so-called 
Mandelstam relations [34]), and the Jacobian is huge. The reformulation 
was done [35] in the language of Schwinger-Dyson or loop equations which 
will be described in the next Section. 

Schwinger-Dyson equations are a convenient way of performing the 
semiclassical expansion, which is an alternative to the path integral. Let 
us illustrate an idea how to do this by an example of the theory. The 
RHS of the Schwinger-Dyson equations is proportional to the Planck's con- 
stant U. In the semiclassical limit ^ — 0, we get 

(-dl + m^) {<p {Xi) ...if {xn)) + ^ (<^2 (xi) ...if {xn)) = , (121) 
whose solution is of the factorized form 

{ip (Xi) ...if [Xn)) = {if (xi)) ...{if {Xn)) + O {h) (122) 

provided that 



obeys 



{<f{x)) = ^,,{x) (123) 

(_a2 + ^2 j ^ A ^2^ ^ Q _ ^^24) 

Equation (124) is nothing but the classical equation of motion for the 
(p^ theory, which specifies extrcma of the action entering the path integral. 
Thus, we have reproduced, using the Schwinger-Dyson equations, the well- 
known fact that the path integral is dominated by a classical solution as 
7i ^ 0. It is also clear how to perform the semiclassical expansion in fi in 
the language of the Schwinger-Dyson equations: one should solve them by 
iterations. 

The reformulation of multicolor QCD in terms of the loop functionals 
$ (C) is, in a sense, a realization of the idea of the master field in the weak 
sense, when the master field acts as an operator in the space of loops. 
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Remark on the large-Nc limit as statistical averaging 

There is yet another, pure statistical, explanation why the large-iVc limit is a 

"scmiclassical" limit for the collective variables $ (C). The matrix V^^ [Cxx], 
that describes the parallel transport along a closed contour Cxxi can be 
reduced by the gauge transformation to 



U [Cxo:] = n [Cxx] diag ( e ^^"^ ^^\...,e ^^"^c (C)^ [Cxx] ■ (125) 
Then $ (C) reads 



The phases aj{C) are gauge invariant and normalized so that aj{C) ~ 1 
as Nc ^ oo. For simplicity we omit below all the indices (including space 
ones) except color. 

The commutator of <I>'s can be estimated using the representation (126). 
Since [a^, aj] oc 6ij , one gets 



in the limit (71), i.e. the commutator can be neglected as Nc — ^ oo, and 
the field $ (C) becomes classical. 

Note that the commutator (127) is of order l/N^- One factor 1/Nc is 
because of g in the definition (126) of <I> (C), while the other has a deep 
reason. Let us image the summation over j in Eq. (126) as some statistical 
averaging. It is well-known in statistics that such averages weakly fluctuate 
as Nc — > oo, so that the dispersion is of order l/Nc- It is the factor which 
emerges in the commutator (127). 

We see that the factorization is valid only for the gauge-invariant quanti- 
ties which involve the averaging over the color indices, like that in Eq. (126). 
There is no reason to expect factorization for gauge invariants which do not 
involve this averaging, for instance for the phases aj{C). Moreover, their 
commutator is ~ 1, so that aj(C)'s strongly fluctuate even &t Nc = oo. 
An explicit example of such strongly fluctuating gauge- invariant quantities 
was first constructed in Ref. [30]. 

The resume to this Remark is that the factorization is due to the ad- 
ditional statistical averaging in the large- A'^c limit. There is no reason to 
assume an existence of the master field in the strong sense in order to 
explain the factorization. 




(126) 



[$(C),<I>(CO] 



1 1 



(127) 



Nc N^ 
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4. QCD in Loop Space 

QCD can be entirely reformulated in terms of the colorless composite field 
$ (C) — the trace of the Wilson loop for closed contours. This fact involves 
two main steps: 

i) All the observables are expressed via $ (C). 

ii) Dynamics is entirely reformulated in terms of $ (C). 

This approach is especially useful in the large- A'c limit where everything 
is expressed via the vacuum expectation value of # (C) — the Wilson loop 
average. Observables are given by summing the Wilson loop average over 
paths with the same weight as in free theory. The Wilson loop average 
obeys itself a close functional equation — the loop equation. 

We begin this Section with presenting the formulas which relate observ- 
ables to the Wilson loops. Then we translate quantum equation of motion 
of Yang-Mills theory into loop space. We derive the closed equation for the 
Wilson loop average as A^^ — ^ oo and discuss its various properties, includ- 
ing a non-perturbative regularization. Finally, we briefly comment on what 
is known about solutions of the loop equation. 

4.1. OBSERVABLES IN TERMS OF WILSON LOOPS 

All observables in QCD can be expressed via the Wilson loops ^> (C) defined 
by Eq. (120). This property was first advocated by Wilson [36] on a lattice. 
Calculation of QCD observables can be divided in two steps: 

1) Calculation of the Wilson loop averages for arbitrary contours. 

2) Summation of the Wilson loop averages over the contours with some 
weight depending on a given observable. 

At finite Nc, observables are expressed via the n-loop averages 



which are analogous to the n-point Green functions for tp^ theory. The 

appropriate formulas for the continuum theory can be found in Ref. [37]. 

Great simplifications occur in these formulas at A^c = oo, when all ob- 
servables are expressed only via the one-loop average 



This is associated with the quenched approximation. 

For example, the average of the product of two colorless quark vector 
currents (101) is given at large Nc by 



(128) 




(129) 



(130) 



C^X1,X2 
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a) b) 

Figure 14- Contours in the sum over paths representing observables: a) in Eq. (130) 
and b) in Eq. (131). The contour a) passes two nailed points xi and X2. The contour b) 
passes three nailed points xi, X2, and xs. 



where the sum runs over contours C passing through the points xi and X2 
as is depicted in Figure 14a. An analogous formula for the (connected) 
correlators of three quark scalar currents reads 

(V^V {X2) V^V (^3)),o,n = E J{C){^{C)) , (131) 

C3X1,X2,XZ 

where the sum runs over contours C passing through the three points xi, 
X2, and X3 as is depicted in Figure 14b. A general (connected) correlator 
of n quark currents is given by a similar formula with C passing through n 
points Xl, . . . , Xn (some of them may coincide). 

The weights J^,, (C) in Eq. (130) and J (C) in Eq. (131) are completely 
determined by free theory. If quarks were scalars rather than spinors, then 
we would get 



J (C) = e -^""'^-Uo <^*^2W = e -'"•^(^) I scalar quarks I , (132) 



where L(C) is the length of the (closed) contour C. For spinor quarks, an ad- 
ditional disentangling of the gamma-matrices is needed (ee, e.g., Ref. [38]). 

Remark on renormalization of Wilson loops 

Perturbation theory for W{C) can be obtained by expanding the path- 
ordered exponential in the definition (129) in g (sec Eq. (84)) and averaging 
over the gluon field A^. Because of ultraviolet divergencies, we need a (gauge 
invariant) regularization. After such a regularization introduced, the Wilson 
loop average for a smooth contour C of the type in Figure 15a reads 

2 (JV^-l) L(C) 

W{C)=e — W^ren (C) , (133) 

where a is the cutoff, L{C) is the length of C, and Wren {C) is finite when 
expressed via the renormalized charge The exponential factor is due to 
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a) b) 

Figure 15. Examples of a) smooth contour and b) contour with a cusp. The tangent 
vector to the contour jumps through angle 7 at the cusp. 

the renormalization of the mass of a heavy test quark. This factor does not 
emerge in the dimensional regularization where d = A—e. The multiphcative 
renormahzation of the smooth Wilson loop was shown in Refs. [39, 40, 41]. 

If the contour C has a cusp (or cusps) but no self-intersections as is 
illustrated by Figure 15b, then W {C) is still multiplicatively renormaliz- 
able [42]: 

VF(C) = Z(7) I^ren(C), (134) 

while the (divergent) factor Z (7) depends on the cusp angle (or angles) 7 
(or 7's) and VFren {C) is finite when expressed via the renormalized charge 

4.2. SCHWINGER-DYSON EQUATIONS FOR WILSON LOOP 

Dynamics of (quantum) Yang-Mills theory is described by the quantum 
equation of motion 

which is understood in the weak sense, i.e. for the averages 

{ ^fFl (x) Q[A]) = n{j^^Q [A]j . (136) 

The standard set of Schwinger-Dyson equations of Yang-Mills theory 
emerges when the functional Q [A] is chosen in the form of the product 
of ^'s as in Eq. (82). 

Strictly speaking, the last statement is incorrect, since we have not 
added, in Eqs. (135) and (136), contributions coming from the variation 
of gauge-fixing and ghost terms in the Yang-Mills action. However, these 
two contributions are mutually cancelled for gauge-invariant functionals 
Q[A]. We shall deal below only with such gauge- invariant functionals (the 
Wilson loops). This is why we have not considered the contribution of the 
gauge-fixing and ghost terms. 
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Figure 16. Contours Cyx and Cxy which enter the RHS's of Eqs. (137) and (141). 

It is also convenient to use the matrix notation (57), when Eq. (135) for 
the Wilson loop takes on the form 

(137) 

where we have restored the units with h = 1. 

The variational derivative on the RHS can be calculated by virtue of 
the formula 

XAij (y) / 1 \ 

= V 6('^ {X - y) [8H>^^ - -Sn^^) (138) 
which is a consequence of 



6AI {x) 



The second term in the parentheses in Eq. (138) — same as in Eq. (64) — 
is because is a matrix from the adjoint representation of SU(A'c). 

By using Eq. (138), we get for the variational derivative on RHS of 
Eq. (137): 



tr— ^Pe^^c<^S^^M = iidy,5^'^Ux-y) 

OAy (X) J 



c 



— trPe — trPe '"^f ^trPe 

Nr Nr 



■ (140) 



The contours Cyx and Cxy, which are depicted in Figure 16, are the parts 
of the loop C: from a; to y and from y to x, respectively. They are always 
closed due to the presence of the delta-function. It implies that x and y 
should be the same points of space but not necessarily of the contour {i.e. 
they may be associated with different values of the parameter a). 
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Wc finally rewrite Eq. (137) as 



Nc 



c 

where we have introduced 



(141) 



A = ^. (142) 

Notice that the RHS of Eq. (141) is completely represented via the (closed) 
Wilson loops. 

4.3. PATH AND AREA DERIVATIVES 

As wc already mentioned, the RHS of Eq. (141) is completely represented 
via the (closed) Wilson loops. It is crucial for the loop-space formulation of 
QCD that the LHS of Eq. (141) can also be represented in loop space as 
some operator applied to the Wilson loop. To do this we need to develop a 
differential calculus in loop space. 

Loop space consists of arbitrary continuous closed loops, C. They can 
be described in a parametric form by the functions x^(c7) G L2,^ where 
ctq < (7 < af and /i = 1, . . . , d, which take on values in a c/-dimensional Eu- 
clidean space. The functions x^{o) can be discontinuous, generally speak- 
ing, for an arbitrary choice of the parameter a. The continuity of the loop 
C implies a continuous dependence on parameters of the type of proper 
length. 

The functions x^{a) G L2 which are associated with the elements of 
loop space obey the following restrictions: 

i) The points u = and u = Uf are identified: x^((To) = Xf^{af) — the 
loops are closed. 

ii) The functions x^{a) and A^i,Xi^{a) +a^, with A^i, and independent 
of a, represent the same element of the loop space — rotational and 
translational invariance. 

iii) The functions x^{a) and x^{a') with a' = f{cr), f'{cr) > describe 
the same loop — reparametrization invariance. 

An example of functionals which are defined on the elements of loop space is 
the Wilson loop average (129) or, more generally, the n-loop average (128). 

The differential calculus in loop space is built out of the path and area 
derivatives. 

^Let us remind that I/2 stands for the Hilbert space of functions Xij,{a) whose square 
is integrable over the Lebesgue measure: JJ^ daxf^ia) < 00. 
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The area derivative of a functional (C) is defined by the difference 



1 




(143) 

where an infinitesimal loop SCni,{x) is attached to a given loop at the point 



X in the /iz^-plane and {Sa^i^l stands for the area enclosed by the 5Cf 



For a rectangular loop 5Cij_i,{x), one gets 5a = dx^ A dxn, where the 
symbol A implies antisymmetrization. 

Analogously, the path derivative is defined by 




where 5x^ is an infinitesimal path along which the point x is shifted from 
the loop and \5Xfj\ stands for the length of the dx^. 

These two differential operations are well-defined for so-called function- 
al of the Stokes type which satisfy the backtracking condition they do 
not change when an appendix passing back and forth is added to the loop 
at some point x: 




(145) 



This condition is equivalent to the Bianchi identity of Yang-Mills theory 
and is obviously satisfied by the Wilson loop (129) due to the properties of 
the non-Abelian phase factor. Such functionals are known in mathematics 
as Chen integrals. 

A simple example of the Stokes functional is the area of the minimal 
surface, A^i^i{C). It obviously satisfies Eq. (145). Otherwise, the length 
L{C) of the loop C is not a Stokes functional, since the lengths of contours 
on the LHS and RHS of Eq. (145) arc different. 

For the Stokes functionals, the variation on the RHS of Eq. (143) is 
proportional to the area enclosed by the infinitesimally small loop SC^n, {x) 
and does not depend on its shape. Analogously, the variation on the RHS 
of Eq. (144) is proportional to the length of the infinitesimal path 5Xfj_ and 
does not depend on its shape. 
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If X is a regular point (like any point of the contour for the func- 
tional (129)), the RHS of Eq. (144) vanishes due to the backtracking con- 
dition (145). In order for the result to be nonvanishing, the point x should 
be a marked (or irregular) point. A simple example of the functional with 
a marked point x is 

$«[C,,] = ^tr (^eVe'^c^^'^"^"^^^) (146) 

with the S\J{Nc) generator inserted in the path-ordered product at the 
point X. 

The area derivative of the Wilson loop is given by the Mandelstam 
formula 

^ ^ ■ixVe'^c''^''^'^ = ^i,^F^,{x)e'^c''^'^- . (147) 



5a^^ {x) Nc Nc 

In order to prove it, it is convenient to choose 5Cny{x) to be a rectangle in 
the /xi/-plane and straightforwardly use the definition (143). The sense of 
Eq. (147) is very simple: F^^ is a curvature associated with the connection 
A,. 

The functional on the RHS of Eq. (147) has a marked point x, and is of 
the type in Eq. (146). When the path derivative acts on such a functional 
according to the definition (144), the result reads 



a-i_trP5(x) e'^c'^^'^^ = i_trPV^5(x) e'^c'^^'^^, (148) 
where 



V^B = d^B-i[Ai„B] (149) 

is the covariant derivative in the adjoint representation. 

Combining Eqs. (147) and (148), we finally represent the expression on 
the LHS of Eq. (137) (or Eq. (141)) as 

i.e. via the action of the path and area derivatives on the Wilson loop. It 
is therefore rewritten in loop space. 

A resume of the results of this subsection is presented in Table 2 as 
a vocabulary for translation of Yang-Mills theory from the language of 
ordinary space in the language of loop space. 
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Ordinary space 



Loop space 





phase factor 


<b{C) 


loop functional 


F^. (x) 


field strength 


5 


area derivative 




covariant derivative 




path derivative 


VAF = 


Bianchi identity 




Stokes functionals 


= 5/5A^ 


Schwinger-Dyson 
equations 




Loop 
equations 



TABLE 2. Vocabulary for translation of Yang-Mills theory from 
ordinary space in loop space. 



Remark on Bianchi identity for Stokes functionals 

The backtracking relation (145) can be equivalently represented as 

by choosing the appendix in Eq. (145) to be an infinitesimal straight line 
in the p-direction and geometrically applying the Stokes theorem. Using 
Eqs. (147) and (148), Eq. (151) can in turn be rewritten as 




(152) 



Therefore, Eq. (151) represents the Bianchi identity in loop space. 

Remark on relation to variational derivative 

The standard variational derivative, 6/6xi^{a), can be expressed via the 
path and area derivatives by the formula 

sr cm 
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where the sum on the RHS is present for the case of a functional having m 
marked (irregular) points Xi = x{ai). A simplest example of the functional 
with m marked points is just a function of m variables xi, . . . , x^- 

By using Eq. (153), the path derivative can be calculated as the limiting 
procedure 

O-+0 



The result is obviously nonvanishing only when is applied to a functional 
with x{a) being a marked point. 

It is nontrivial that the area derivative can also be expressed via the 
variational derivative [40]: 



(T+O 

^ / da'(a' -a)—^—^. (155) 



(5(T^^(x((7)) 

(T — 

The point is that the six-component quantity, 6/da^i^{x{a)), is expressed 
via the four-component one, S/dXfj,{a), which is possible because the com- 
ponents of 5/5ani,{x{a)) are dependent due to the loop-space Bianchi iden- 
tity (151). 

4.4. LOOP EQUATIONS 

By virtue of Eq. (150), Eq. (141) can be represented completely in loop 
space: 



$(C) 



c 

or, using the definitions (128) and (129) of the loop averages, as 



(156) 



W{C) = xj dy, 5^"^ {x - y) [w^ {Cy,,C,y) - {C) 
c " 



(157) 

This equation is not closed. Having started from W{C), wc obtain an- 
other quantity, VF2 (C*!, (^2), so that Eq. (157) connects the one-loop av- 
erage with a two-loop one. This is similar to the case of the (quantum) 
(^^-theory, whose Schwinger-Dyson equations connect the n-point Green 
functions with different n. We shall derive this complete set of equations 
for the n-loop averages in this Subsection later on. 



LARGE-N GAUGE THEORIES 



47 



However, the two-loop average factorizes in the large- A^c limit: 

W2{Ci,C2) = W{Ci)W{C2)+o(^^y (158) 

as was discussed in Subsection 3.4. Keeping the constant A (defined by 
Eq. (142)) fixed in the large- AT^ limit as is prescribed by Eq. (71), we get [35] 

^^ 6a l(x) ^ = ^fdy, S^''^ ix-y)W (Cy,) W (C^y) (159) 
as iVc — *• oo. 

Equation (159) is a closed equation for the Wilson loop average in the 
large- iVc limit. It is referred to as the loop equation. 

To find W{C), Eq. (159) should be solved in the class of Stokes func- 
tional with the initial condition 

W{0) = 1 (160) 

for loops which are shrunk to points. 

The factorization (158) can itself be derived from the chain of loop 
equations. Proceeding as before, we get 

A ^d(7^„{x) 

1 



^n+l {Cxyj Cyx, . . . , Cji) — 'j;^^n (Clj • • • ) C'n 



= jdy^S^'^Hx-y) 

Ci 



j>2 "'c 



Wn{Ci,...,Cn) ■ (161) 



Here x belongs to Ci; CiCj stands for the joining of Ci and Cj; Cj means 
that Cj is omitted. 

Equation (161) looks like the Schwinger-Dyson equation for the tp^- 
theory. Moreover, the number of colors Nc enters Eq. (161) simply as a 
scalar factor iV~^, likewise Plank's constant h enters in the -theory. It is 
the major advantage of the use of loop space. What is said in Subsection 3.6 
about the "semiclassical" nature of the 1/iVc-expansion of QCD is explicitly 
realized in Eq. (161). Its expansion in 1/Nc is straightforward. 

At A^c = CO, Eq. (161) is simplified to 

^ -Wn{Ci,...,Cn) = \j> dy, 5('^)(x - y)Wn+l{Cyx, Cxy, . . . , Cn). 



5cT^^ (x) 

Ci 



(162) 
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Figure 17. Graphic representation of the terms on the RHS of Eq. (164). 
This equation possesses a factorized solution 



provided W (C) obeys Eq. (159) which plays the role of a "classical" equa- 
tion in the large- A^c limit. Thus, we have given a non-perturbative proof of 
the large-iVc factorization of the Wilson loops. 

4.5. RELATION TO PLANAR DIAGRAMS 

The perturbation-theory expansion of the Wilson loop average can be cal- 
culated from Eq. (84) which we represent in the form 



where 9cil,2, . . . ,n) orders the points xi, . . ., Xn along contour in the cyclic 
order and is given by Eq. (82). This ^-function has the meaning of 

the propagator of a test heavy particle which lives in the contour C. 

We assume, for definitivcness, the dimensional regularization through- 
out this Subsection to make all the integrals well-defined. 

Each term on the RHS of Eq. (164) can be conveniently represented 
by the diagram in Figure 17, where the integration over the contour C is 
associated with each point Xi lying in the contour C. 

These diagrams are analogous to those discussed in Subsection 3.2 with 
one external boundary — the Wilson loop in the given case. In the large- A^c 



W„{Cu...,C„) = W{C,)---W(Cn)+0 




(164) 
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a) b) c) 

Figure 18. Planar diagrams for W(C): a) of order A with gluon propagator, and of order 
b) with two noninteracting gluons and c) with the three-gluon vertex. Diagrams of 
order A^ with one-loop insertions to gluon propagator are not drawn. 

limit, only planar diagrams survive. Some of them, which are of the lowest 
order in A, are depicted in Figure 18. 

The large- iVc loop equation (159) describes the sum of the planar dia- 
grams. Its iterative solution in A reproduces the set of planar diagrams for 
W (C) provided the initial condition (160) and some boundary conditions 
for asymptotically large contours are imposed. 

Equation (164) can be viewed as an ansatz for W (C) with some un- 
known functions G^"---//„ (a:;i, • • • ,Xn) to be determined by the substitution 
into the loop equation. To preserve symmetry properties of W (C), the 
functions G^^^ must be symmetric under a cyclic permutation of the points 
1, . . ., n and depend only on Xi — Xj (translational invariancc). A main 
advantage of this ansatz is that it automatically corresponds to a Stokes 
functional, due to the properties of vector integrals, and the initial condi- 
tion (160) is satisfied. 

The action of the area and path derivatives on the ansatz (164) is easily 
calculable. For instance, the area derivative reads 



5W{C) 



= £ /dxf ... i'd<"^e(l,2,...,n) 

c c 



(n+l) 



P-n (^) Xl, . . . , Xji) 



+ {SnpS„a - G^^p^l,,^^ {z,z,xi,..., x^ 



(165) 



The analogy with the Mandelstam formula (147) is obvious. 

More about solving the loop equation by the ansatz (164) can be found 
in Refs. [37, 43, 44]. 
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4.6. LOOP-SPACE LAPLACIAN AND REGULARIZATION 

The loop equation (159) is not yet entirely formulated in loop space. It is a 
d-vector equation whose both sides depend explicitly on the point x which 
does not belong to loop space. The fact that wc have a d-vector equation 
for a scalar quantity means, in particular, that Eq. (159) is overspecified. 

A practical difficulty in solving Eq. (159) is that the area and path 
derivatives, 5/5a^j^v{x) and 5^, which enter the LHS are complicated, gen- 
erally speaking, non-commutative operators. They are intimately related 
to the Yang-Mills perturbation theory where they correspond to the non- 
Abelian field strength F^i, and the covariant derivative V^. However, it is 
not easy to apply these operators to a generic functional W{C) which is 
defined on elements of loop space. 

A much more convenient form of the loop equation can be obtained 
by integrating both sides of Eq. (159) over dxn along the same contour C, 
which yields 



j> dxv d'^ 



^ ^(^) = ^jdx^j dy^ 5^''^ {x - y) W{Cy,)W{C,y) 



c c c 

(166) 

Now both the operator on the LHS and the functional on the RHS are 
scalars without labeled points and are well-defined in loop space. The oper- 
ator on the LHS of Eq. (166) can be interpreted as an infinitesimal variation 
of elements of loop space. 

Equations (159) and (166) are completely equivalent. A proof of equiv- 
alence of scalar Eq. (166) and original d-vector Eq. (159) is based on the 
important property of Eq. (159) whose both sides are identically annihilated 
by the operator df, . It is a consequence of the identity 

V^V. F^, = F^,] = (167) 

in the ordinary space. Due to this property, the vanishing of the contour 
integral of some vector is equivalent to vanishing of the vector itself, so that 
Eq. (159) can in turn be deduced from Eq. (166). 

Equation (166) is associated with the so-called second-order Schwinger- 
Dyson equation 

/ Av,f;„ (.) ^ - nf Myii'^ (x - ,) 

(168) 

in the same sense as Eq. (159) is associated with Eq. (135). It is called "sec- 
ond order" since the RHS involves two variational derivatives with respect 
to Aj,. 
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The operator on the LHS of Eq. (166) is a well-defined object in loop 
space. When applied to regular functional which do not have marked 
points, it can be represented, using Eqs. (154) and (155), in an equiva- 
lent form 



C 



da / da' 



(TO a-O 



6x^{a') dxi^{a) 



(169) 



As was first pointed out by Gervais and Neveu [45], this operator is noth- 
ing but a functional extension of the Laplace operator, which is known in 
mathematics as the Levy operator.^ Equation (166) can be represented in 
turn as an (inhomogeneous) functional Laplace equation 

AW{C) = X^dx^f dy^ S^'^^ {x - y) W{Cy^)W{C^y) . (170) 



We shall refer to this equation as the loop-space Laplace equation. 

The form (170) of the loop equation is convenient for a non-perturbative 
ultraviolet regularization. 

The idea is to start from the regularized version of Eq. (168), replacing 
the delta-function on the RHS by the kernel of the regularizing operator: 



5-^5^''\x-y)^(y 



R'^'^ 5^'^\x - y) 



with 



R 



ab 



ab 



(171) 



(172) 



where is the covariant derivative in the adjoint representation. The 
regularized version of Eq. (168) is 



y \y 



R 



ab 



SAb (y) 6 A- (x) ■ 

(173) 

To translate Eq. (173) in loop space, we use the path-integral represen- 
tation 



R 



ab 



2 

= I Dr{t) e-'^ lo "''"^'^ 2tT: [eU{ry,)t'U{r,y)] (174) 



r(0) = x 



with 



C/(r 



yx) 



Pe 



(175) 



^See the book by Levy [46] and a review [47]. 



52 



YURI MAKEENKO 




Figure 19. Contours Cy^Vxy and CxyVyx which enter the RHS's of Eqs. (176) and (177). 

where the integration is over regulator paths from x to y whose typical 
length is ~ a. 

Calculating the variational derivatives on the RHS of Eq. (173), using 
Eq. (174) and the completeness condition (64), we get as N ^ oo: 

c c 

2 

X j Vrit) e-Uo W $ {Cy,r,y) $ (C,^r^.) , (176) 

r(0) = x 

where the contours Cyxrxy and CxyVyx arc depicted in Figure 19. Averaging 
over the gauge field and using the large- A'^c factorization, we arrive at the 
regularized loop-space Laplace equation [48] 

2 

AW{C) = X^dx^^ dy^ j Dr{t) e io" W{Cyxrxy)W{Cxyryx) 

C C r(0) = x 

r(a^) — y 

(177) 

which manifestly recovers Eq. (170) when a ^ 0. 

The constructed regularization is non-perturbative while perturbatively 
reproduces regularized Feynman diagrams. An advantage of this regular- 
ization of the loop equation is that the contours CyxVxy and CxyVyx on the 
RHS of Eq. (177) both are closed and do not have marked points if C does 
not have. Therefore, Eq. (177) is written entirely in loop space. 

Remark on functional Laplacian 

It is worth noting that the representation of the functional Laplacian on 
the RHS of Eq. (169), which involves the standard variational derivatives, 
is defined for a wider class of functionals than Stokes functionals. It is easier 
to deal with the whole operator A, rather than separately with the area 
and path derivatives. 
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The functional Laplacian is parametric invariant and possesses a num- 
ber of remarkable properties. While a finite-dimensional Laplacian is an 
operator of the second order, the functional Laplacian is that of the first 
order and satisfies the Leibnitz rule 

A (UV) = A{U)V + UA {V) . (178) 

The functional Laplacian can be approximated [49] in loop space by a 
(second-order) partial differential operator in such a way to preserve these 
properties in the continuum limit. This loop-space Laplacian can be in- 
verted to determine a Green function G (C, C") in the form of a sum over 
surfaces Sc,c' connecting two loops which is analogous to the sum-over- 
path representation of the Green function of the ordinary Laplacian. The 
standard perturbation theory can then be recovered by iterating Eq. (170) 
(or its regularized version (177)) in A with the Green function of the loop- 
space Laplacian. 

4.7. SURVEY OF NON-PERTURBATIVE SOLUTIONS 

While the loop equations were proposed long ago, not much is known about 
their non-perturbativc solutions. We briefly list some of the results. 
It was shown in Ref. [50] that area law 

W (C) = (C)) cx e -K-A^UC) (179) 

satisfies the large- A/'c loop equation for asymptotically large C. However, a 
self-consistency equation for K, which should relate it to the bare charge 
and the cutoff, was not investigated. In order to do this, one needs more 
detailed information about the behavior of W (C) for intermediate loops. 

The free bosonic Nambu-Goto string which is defined as a sum over 
surfaces spanned by C 

W{C)= e-^-^(^), (180) 

S:dS=C 

with the action being the area A (S) of the surface S, is not a solution 
for intermediate loops. Consequently, QCD docs not reduce to this kind 
of string, as was originally expected in Refs. [51, 52, 53, 54, 55]. Roughly 
speaking, the ansatz (180) is not consistent with the factorized structure 
on the RHS of Eq. (159). 

Nevertheless, it was shown that if a free string satisfies Eq. (159), then 
the same interacting string satisfies the loop equations for finite Nc- Here 
"free string" means, as usual in string theory, that only surfaces of genus 
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zero are present in the sum over surfaces, while surfaces or higher gen- 
era are associated with a string interaction. The couphng constant of this 
interaction is O . 

A formal solution of Eq. (159) for all loops was found by Migdal [56] in 
the form of a fcrmionic string 

W{C)= J2 /"L>V'e-/'^'^[^^^^'=^+^'^'"'^], (181) 

S:dS=C 

where the world sheet of the string is parametrized by the coordinates 
^1 and ^2 for which the 2-dimensional metric is conformal, i.e. diagonal. 
The field V'(0 describes 2-dimensional elementary fermions (elves) living 
in the surface S, and m stands for their mass. Elves were introduced to 
provide factorization which now holds due to some remarkable properties 
of 2-dimcnsional fermions. For large loops, the internal fermionic structure 
becomes frozen, so that the empty string behavior (179) is recovered. For 
small loops, the elves are necessary for asymptotic freedom. However, it 
is unclear whether or not the string solution (181) is practically useful for 
a study of multicolor QCD, since the methods of dealing with the string 
theory in four dimensions are not yet developed. 

A very interesting solution of the large-iVc loop equation on a lattice 
was found by Eguchi and Kawai [57]. They showed that the SU(A'c) gauge 
theory on an infinite lattice reduces at Nc = oo to the model on a hyper- 
cube. The equivalence is possible only at A^c = oo, when the space-time 
dependence is absorbed by the internal symmetry group. More about this 
large- ATp reduction will be said in the next Section. 

4.8. WILSON LOOPS IN QCD2 

Two-dimensional QCD is popular since the paper by 't Hooft [58] as a 
simplified model of QCD4. 

One can always choose the axial gauge Ai = 0,so that the commutator in 
the non-Abelian field strength (60) vanishes in two dimensions. Therefore, 
there is no gluon self-interaction in this gauge and the theory looks, at the 
first glance, like the Abelian one. 

The Wilson loop average in QCD2 can be straightforwardly calculated 
via the expansion (164) where only disconnected (free) parts of the corre- 
lators G^") for even n should be left, since there is no interaction. Only the 
planar structure of color indices contributes at AT^ = 00. Diagrammatically, 
the diagrams of the type depicted in Figure 18a and Figure 18b are relevant 
for contours without self-intersections, while that in Figure 18c should be 
omitted in two dimensions. 



LARGE-N GAUGE THEORIES 



55 



X2 = 2/2l 




Figure 20. Graphic representation of the contour integral on the LHS of Eq. (186) in 
the axial gauge. The bold line represents the gluon propagator (184) with X2 = y2 due 
to the delta-function. 



The color structure of the relevant planar diagrams can be reduced by 
using the completeness condition (64) at large Nc- We have 




x6lc(l, 2,..., 2k) Df,^y^ {xi - X2) • • • D^_,^y^ {x2k-i ~ X2k) , (182) 

where the points xi, . . ., X2k are still cyclic ordered along the contour. We 
can exponentiate the RHS of Eq. (182) to get finally 

W{C)= e '^^'^ §c dy^D^^ix-y) _ ^^g3^| 

This is the same formula as in the Abelian case if A stands for e^. 

The propagator D^^^ (x, y) is, strictly speaking, the one in the gauge 
^1=0 which reads 

D^^ {x-y) = ^d^2Su2 \xi - yi\ 6^^^ (x2 - y2) ■ (184) 

However, the contour integral on the RHS of Eq. (183) is gauge invariant, 
and we can simply choose instead 

1 ^2 

D^i-u {x-y) = 5u,y — In , (185) 

47r {x - y) 

where I is an arbitrary parameter of the dimension of length. Nothing 
depends on it because the contour integral of a constant vanishes. 

The contour integral in the exponent on the RHS of Eq. (183) can be 
graphically represented as is depicted in Figure 20, where X2 = 2/2 due to 
the delta-function in Eq. (184) and the bold line represents \x\ — yi\. This 
gives 

jdx^" jdy''D^^{x-y) = A{C) (186) 
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a) 



b) 



Figure 21. Contours with one self- intersection: Ai and A2 stand for the areas of the 
proper windows. The total area enclosed by the contour in Figure a) is Ai-\-A2. The areas 
enclosed by the exterior and interior loops in Figure b) are A\ + A2 and A2, respectively, 
while the total area of the surface with the folding is Ai + 2A2. 

where A (C) is the area enclosed by the contour C. We get finally 



for the contours without self-intersections. 

Therefore, area law holds in two dimensions both in the non-Abelian 
and Abelian cases. This is, roughly speaking, because of the form of the 
two-dimensional propagator (185) which falls down with the distance only 
logarithmically in the Feynman gauge. 

The difference between the Abelian and non-Abelian cases shows up for 
the contours with self-intersections. 

We first note that the simple formula (186) does not hold for contours 
with arbitrary self-intersections. 

The simplest contours with one self-intersection are depicted in Fig- 
ure 21. There is nothing special about the contour in Figure 21a. Equa- 
tion (186) still holds in this case with A (C) being the total area A (C) = 
A1 + A2. 

The Wilson loop average for the contour in Figure 21a coincides both 
for the Abelian and non-Abelian cases and equals 



W{C)= e-t^(^) 



(187) 



W{C) 



= e 



t(Al+A2) 



(188) 



This is nothing but the exponential of the total area. 
For the contour in Figure 21b, we get 




(189) 



c 



c 



This is easy to understand in the axial gauge where the ends of the propa- 
gator line can lie both on the exterior and interior loops, or one end at the 
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a) b) c) d) 

Figure 22. Three type of contribution in Eq. (189) The ends of the propagator Une 
lie both on a) exterior and b) interior loops, or c), d) one end on the exterior loop and 
another end on the interior loop. 

exterior loop and the other end on the interior loop. These cases are illus- 
trated by Figure 22. The contributions of the diagrams in Figure 22a,b,c,d 
are A\ + A2, A2, A2, and A2, respectively. The result given by Eq. (189) is 
obtained by summing over all four diagrams. 

For the contour in Figure 21b, the Wilson loop average is 

11/(C) = 6-5(^1+4^2) (190) 

in the Abelian case and 

VF(C) = (1- Ayla) e-t(^i+2^2) (191) 

in the non- Abelian case at A^c = 00. They coincide only to the order A 
as they should. The difference to the next orders is because only the dia- 
grams with one propagator line connecting the interior and exterior loops 
are planar and, therefore, contribute in the non-Abelian case. Otherwise, 
the diagram is non-planar and vanishes as AT^ — 00. Notice, that the expo- 
nential of the total area A (C) = ^1 + 2A2 of the surface with the folding, 
which is enclosed by the contour C, appears in the exponent for the non- 
Abelian case. The additional pre-exponential factor could be associated 
with an entropy of foldings of the surface. 

The Wilson loop averages (188) and (191) in QCD2 at large Nc as well as 
the ones for contours with arbitrary self-intersections, which have a generic 
form 

W(C) = P(Ai,...,^„)e-i^^^« (192) 

where P is a polynomial of the areas of individual windows and Area is 
the total area of the surface with foldings, were first calculated in Ref. [59] 
by solving the two-dimensional loop equation and in Ref. [60] by applying 
the non-Abelian Stokes' theorem. 
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Remark on the string representation 

A nice property of QCD2 at large is that the exponential of the area en- 
closed by the contour C emerges'^ for the Wilson loop average W (C). This 
is as it should for the Nambu-Goto string (180). However, the additional 
pre-exponential factors (like that in Eq. (191)) are very difficult to interpret 
in the stringy language. They may become negative for large loops which 
is impossible for a bosonic string. This explicitly demonstrates in d = 2 the 
statement of the previous subsection that the Nambu-Goto string is not a 
solution of the large- A'^c loop equation. 

5. Large-AT Reduction 

The large- ATc reduction was first discovered by Eguchi and Kawai [57] who 
showed that the Wilson lattice gauge theory on a d-dimensional hypercu- 
bic lattice is equivalent at = 00 to the one on a hypcrcube with peri- 
odic boundary conditions. This construction is based on an extra (Z^^)'^- 
symmetry which the reduced theory possesses to each order of the strong 
coupling expansion. 

Soon after it was recognized that a phase transition occurs in the re- 
duced model with decreasing the coupling constant, so that this symme- 
try is broken in the weak coupling regime. To cure the construction at 
weak coupling, the quenching prescription was proposed by Bhanot, Heller 
and Neuberger [61] and elaborated by many authors. An elegant alterna- 
tive reduction procedure based on twisting prescription was advocated by 
Gonzalez- Arroyo and Okawa [62] . Each of these prescriptions results in the 
reduced model which is fully equivalent to multicolor QCD, both on the 
lattice and in the continuum. 

While the reduced models look as a great simplification, since the space- 
time is reduced to a point, they still involve an integration over d infinite 
matrices which is in fact a continual path integral. It is not clear at the 
moment whether or not this is a real simplification of the original theory 
which can make it solvable. Nevertheless, the reduced models are useful 
and elegant representations of the original theory at large Nc- 

We shall start this Section by a simplest example of a pure matrix 
scalar theory. The quenched reduced model for this case was proposed by 
Parisi [63] on the lattice end elaborated by Gross and Kitazawa [64] in the 
continuum, while the twisted reduced model was advocated by Eguchi and 
Nakayama [65]. Then we concentrate on the Eguchi-Kawai reduction of 
Yang-Mills theory. 

''This is not true, as is already discussed, in the Abehan case for contours with self- 
intersections. 
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5.1. REDUCTION OF SCALAR FIELD 

Let us begin with a simplest example of a pure matrix scalar theory on a 
lattice whose partition function is defined by the path integral 

^ = / n n d'P^ e ("'''''^'^^'' ''^^^""') . (193) 

a; i>j 

Here (fx is a A^c x -^c Hermitean matrix field with x running over sites of a 
hypercubic lattice and V[(p] is some interaction potential, say 

V[^] = + + (194) 

The prescription of the large- A^c reduction is formulated as follows. We 
substitute 

ip^ ^ S^^Sl, (195) 

where 

[Sx]^i = e ^Pfe^M^/^i = diag ( e 'p^^""^ , . . . , e 'Pn,^'^'^ (195) 

is a diagonal unitary matrix which eats the coordinate dependence, so that 
$ does not depend on x. 

The averaging of a functional F[ipx] which is defined with the same 
weight as in Eq. (193), 



can be calculated at AT^ = 00 by 

"o /i=l 1=1 



- <.« /_: n n lr(^[«-*«ii)a^.„. («8) 



where the average on the RHS is calculated [63] for the quenched reduced 
model whose averages are defined by 



{m 



1 



Reduced ^Reduced 



i>j 

The partition function of the reduced model reads 



^Reduced = / d^., e "^^ ^ ^^^^^^^ ^. l*-!^ ((-^-^^)«) (200) 

i>j 
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which can be deduced, modulo the volume factor, from the partition func- 
tion (193) by the substitution (195). 

Notice that the integration over the momenta on the RHS of Eq. (198) 
is taken after the calculation of averages in the reduced model. Such vari- 
ables are usually called quenched in statistical mechanics which clarifies the 
terminology. 

Since A^c ^ oo it is not necessary to integrate over the quenched mo- 
menta in Eq. (198). The integral should be recovered if 's would be uni- 
formly distributed in a d-dimensional hypercube. Moreover, a similar prop- 
erty holds for the matrix integral over $ as well, which can be substituted 
by its value at the saddle point configuration ^g- 

F[ip,]) ^ F[S,^sSl] , (201) 



where the momenta are uniformly distributed in the hypercube. There- 
fore, this saddle point configuration plays the role of a master field in the 
sense of Subsection 3.5. 

In order to show how Eq. (198) works, let us demonstrate how the 
planar diagrams of perturbation theory for the scalar matrix theory (193) 
are recovered in the quenched reduced model. 

The quenched reduced model (200) is of the general type discussed in 
Section 3. The propagator is given by 

^ij^ky = ^G{pi-pj)6u5kj (202) 



Gauss Nc 



with 



M-E,cos((pf-p^^) ^ 



It is convenient to associate the momenta pi and pj in Eq. (203) with 
each of the two index lines representing the propagator and carrying, re- 
spectively, indices i and j. Remember, that these lines are oriented for a 
Hermitean matrix $ and their orientation can be naturally associated with 
the direction of the flow of the momentum. The total momentum carried 
by the double line is pi — pj. 

The simplest diagram which represents the correction of the second 
order in A3 to the propagator is depicted Figure 23. The momenta pi and 
Pj flows along the index lines i and j while the momentum pk circulates 
along the index line k. The contribution of the diagram in Figure 23 reads 

-§G {pi - pjf Y.G{Pi- Pk) G {pk -Pj), (204) 
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Figure 23. Simplest planar diagram of the second order in A3 for the propagator in 
the quenched reduced model (200). The momentum pi flows along the index line i. The 
momentum — pj is associated with the double line ij. 

where the summation over the index k is just a standard one over indices 
forming a closed loop. 

In order to show that the quenched-model result (204) reproduces the 
correction to the propagator in the original theory on an infinite lattice, we 
pass to the variables of the total momenta flowing along the double lines: 

Pi - Pj =p; Pk-Pj = q; Pi-Pk = p-q, (205) 

which is obviously consistent with the momentum conservation at each 
of the two vertices of the diagram in Figure 23. Since p^'s are uniformly 
distributed in the hypercube, the summation over k can be substituted as 
iVc DO by the integral 

The prescription (198) then gives the correct expression 

«'f^(^)V4(^F^^^)^^^-^) ^'''^ 

for the second-order contribution of the perturbation theory for the prop- 
agator on the lattice. 

It is now clear how a generic planar diagram is recovered by the reduced 
model. We first represent the diagram by the double lines and associate the 
momentum with an index line carrying the index i. Then we write down 
the expression for the diagram in the reduced model with the propaga- 
tor (203). Passing to the momenta flowing along the double lines, similar 
to Eq. (205), we get an expression which coincides with the integrand of the 
Feynman diagram for the theory on the whole lattice. It is crucial that such 
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a change of variables can always be done for a planar diagram consistently 
with the momentum conservation at each vertex. The last step is that the 
summation over indices of closed index lines reproduces the integration over 
momenta associated with each of the loops according to Eq. (206) . It is as- 
sumed that the number of loops is much less than Nc which is always true 
for a given diagram since Nc is infinite. 

We thus have shown how planar diagrams of the lattice theory defined 
by the partition function (193) are recovered by the reduced model (200). 
The lattice was needed only as a regularization to make all integrals well- 
defined and was not crucial in the consideration. This construction can be 
formulated directly for the continuum theory [64, 66] where the propagator 
turns into 

G{Pi-Pj) = (208) 

and a Lorentz- invariant regularization can be achieved by choosing < A^. 

Remark on the twisted reduced model 

An alternative reduction procedure is based on the twisting prescription 
[62]. We again perform the unitary transformation (195) with the matrices 
Sx being expressed via a set of d (unitary) x Nc matrices by 

Sx = r^l/«r^2/aj.^3/aj.^4/a ^209) 

where the coordinates of the (lattice) vector measured in the lattice 

units. The matrices are explicitly constructed in Ref. [62] and commute 

by 

r^r, = v^z^r^ (210) 

with Z^jy = being elements of Ztv^ • 

For the twisting reduction prescription, Eq. (198) is valid providing the 
average on the RHS is calculated for the twisted reduced model which is 
defined by the partition function [65] 



^TRM 



= J d^e-^^^^-^I^l+^^^M f-^''^* . (211) 



We can change the order of E's in Eq. (209) defining a more general 
path-dependent factor 

Sx = P n (212) 

oo 

The path-ordered product in this formula runs over all links I = {z, fi) 
forming a path Cxoo from infinity to the point x. 



LARGE-N GAUGE THEORIES 



63 



Due to Eq. (210), changing the form of the path multiphes Sx by the 
Abehan factor 

Z{C)= n (213) 

neS:dS=C 

where (/x, i^) is the orientation of the plaquette □. The product runs over 
any surface spanned by the closed loop C which is obtained by passing 
the original path forward and the new path backward. Due to the Bianchi 
identity 

n v(°)=i (214) 

□ 6cu&e 

where the product goes over six plaquettes forming a 3-dimensional cube 
on the lattice, the product on the RHS of Eq. (213) does not depend on 
the form of the surface S and is a functional of the loop C. 

It is now easy to see that under this change of the path we get 

[SMSi]ki^\ZiC)\'[SxUSl]M (215) 

and the path-dependence is canceled because |Z(C)|^ = 1. This is a general 
property which holds for the twisting reduction prescription of any even (i.e. 
invariant under the center Zn^) representation of SU{Nc). 

5.2. REDUCTION OF YANG-MILLS FIELD 

The statement of the Eguchi-Kawai reduction of the Yang-Mills field says 
that the theory on a d-dimensional space-time is equivalent at A'^c = oo to 
the reduced model which is nothing but its reduction to a point. The action 
of the reduced model is given by 

SEK = ^^tv[A^,A,]\ (216) 

where ^4^ are d space-independent matrices and A is a dimensionful param- 
eter. 

A naive statement of the Eguchi-Kawai reduction is that the averages 
coincide in both theories, for example, 

/ J_ tr P e ^ / ^^''^-^ A =/±tTFe^f^^'^A (217) 
\Nc /d-dim \Nc /ek ^ ^ 

where the LHS is calculated with the action (59) and the RHS is calculated 
with the reduced action (216). Strictly speaking, this naive statement is 
valid only in d = 2 or supersymmetric case for the reason which will be 
explained in a moment. 
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The precise equivalence is valid only if the average of open Wilson loops 
vanish in the reduced model: 

^t,Pe'-fcy.'^'^^\ =0, (218) 

J^c I EK 

as it does in the d-dimensional theory due to the local gauge invariance 
under which 

/ ij \ / ij 

The point is that this gauge invariance transforms in the reduced model 
into (global) rotation of the reduced field by constant matrices O: 

Af, n'fA^n . (220) 

which does not guarantee such vanishing in the reduced model. 

There exists, however, a symmetry of the reduced action (216) under 
the shift of by a unit matrix^: 

A^^A^ + a^5^^, (221) 

which is often called the R'^-symmetry. Under the transformation (221), we 
get 

(pe'-^'^y^'^'^^y ^ e*(^"--")»'^ (p e'-^'^y^'^'^'^y (222) 

which guarantees, if the symmetry is not broken, the vanishing of the open 
Wilson loops 

WEK{Cy.)^(^tTPe'^^y^'^"^') = (223) 

\ -fVc / EK 

in the reduced model. 

The equivalence of the two theories can then be shown using the loop 
equation which reads for the reduced model 

= (l^trP[A^,[A^,^.]]e^^c..'^«''^'' 
oa^„[x) \l\c / EK 

= XA^ ( — trP— — e ■'Cxx " 

\Nc dA^ /EK 

= XA'^fdy,WEK{Cy^)WEK{C^y). (224) 
c 

*This symmetry is rigorously defined on a lattice where it is associated with a 
direction-dependent ZjVc transformation. 
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The RHS is pretty much similar to the one in Eq. (159) while 6^'^\x - y) 
is missing. 

This delta function can be recovered if the R'^ symmetry is not broken 
since 

Wek (Cy.) ~ -^^^Wek (Cy.) (225) 

due to Eq. (223) for the open loops. 

This is not a rigorous argument since a regularization is needed. What 
actually happens is the following. If we smear the delta function introducing 

S^\x) = (^ye-^'^y\ (226) 



. 2nJ 
then 

-^(<5r(0))^ocA^e-^^^^5W(x), (227) 
reproducing the delta function. 

5.3. R^-SYMMETRY IN PERTURBATION THEORY 

Since Nc is infinite, the R^-symmetry can be broken spontaneously. The 
point is that the large- iVc limit plays the role of a statistical averaging, as 
is mentioned already in Subsection 3.6, and phase transitions are possi- 
ble for infinite number of degrees of freedom. This phenomenon occurs in 
perturbation theory of the reduced model for d > 3. 

The perturbation theory can be constructed expanding the fields around 
solutions of the classical equation 

[Ai„[A^,A,]]=0. (228) 

Any diagonal matrix 

= diag {;>«,..., (229) 

is a solution to Eq. (228). 

The perturbation theory of the reduced model can be constructed ex- 
panding around the classical solution (229): 

A^ = A^i + 5^^, (230) 

where A^ is off-diagonal. 

Substituting (230) into the action (216), we get 

Sek = tr |l[p/., - ^^]'} + higher orders . (231) 
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To fix the gauge symmetry (220), it is convenient to add 

-^g-f. = tr |^bM>^^]' + MW^c]^ , (232) 

where b and c are ghosts. 

The sum of (231) and (232) gives 

S2 = tr + \p^^,b]\p^,c]^ (233) 

up to quadratic order in A'^. 

Doing the Gaussian integral over A^, we get at the one-loop level: 




where the integration over accounts for equivalent classical solutions. 

For d = 1 the product on the RHS of Eq. (234) reproduces the Vander- 
monde determinant. For d = 2 it vanishes and does not affect dynamics. 
For d > 3 the measure is singular and the eigenvalues collapse. This leads 
us to a spontaneous breakdown of the R'^ in perturbation theory. 

The equivalence between the Nc = oo Yang-Mills theory on a whole 
space and the reduced model can be provided [61] introducing a quenching 
prescription similar to the one described in Subsection 5.1. Then no collapse 
of eigenvalues happens and d-dimensional planar graphs are reproduced by 
the rcdiiccd model. More about the quenching prescription in Yang-Mills 
theory can be found in the reviews [44, 67] and cited there original papers. 

Remark on supersymmetric case 

In a supersymmetric gauge theory, there is an extra contribution from 
fermions to the exponent on the RHS of Eq. (234). Since the integration 
over fermions results in the extra factor [(p^*^ — p^^)^] this yields fi- 
nally the exponent l — d/2 + tr 1/2. It vanishes in d = 4 for either Majorana 
or Weyl fermions and in d = 10 for the Majorana-Weyl fermions. There- 
fore, the R'^-symmetry is not broken and no quenching is needed in the 
supersymmetric case [68, 5]. 



5.4. TWISTED REDUCED MODEL 

The continuum version of the twisted reduced model can be constructed [69] 
by substituting Aj^ Aj^ — 7^^^ into the action (216), where the matrices 7^ 
obey the commutation relation 



(235) 
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where B^i, is an antisymmetric tensor and d is even. This is possible only 
for infinite Hermitean matrices (operators). An example of such matrices is 
X and p operators in quantum mechanics. Eq. (235) is a continuum version 
of Eq. (210). 

The Wilson loop averages in the twisted reduced model are defined by 
W^TEK (C,,) = / tr P e '^"^^ 1^ tr P e ' ^^y^ ) . (236) 

\ ^^c J^^c I TEK 

They vanish for open loops which is provided by the vanishing of the trace 
of the path-ordered exponential of in this definition. For closed loops 
this factor does not vanish and is needed to provide the equivalence with 
d-dimensional Yang-Mills perturbation theory, since the classical extrema 
of the twisted reduced model arc A^^ = 7^ and the perturbation theory is 
constructed expanding around this classical solution. 

The proof of the equivalence can be done using the loop equation quite 
similarly to that of Subsection 5.2 for the Eguchi-Kawai model with an 
unbroken R*^ symmetry. 
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